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Abstract. The maximally-decoupled method has been considered as a theory to apply an 
basic idea of an integrability condition to certain multiple parametrized symmetries. The 
method is regarded as a mathematical tool to describe a symmetry of a collective subma- 
nifold in which a canonicity condition makes the collective variables to be an orthogonal 
coordinate-system. For this aim we adopt a concept of curvature unfamiliar in the con- 
ventional time-dependent (TD) self-consistent field (SCF) theory. Our basic idea lies in 
the introduction of a sort of Lagrange manner familiar to fluid dynamics to describe a 
collective coordinate-system. This manner enables us to take a one-form which is linearly 
composed of a TD SCF Hamiltonian and infinitesimal generators induced by collective vari- 
able differentials of a canonical transformation on a group. The integrability condition 
of the system read the curvature C = 0. Our method is constructed manifesting itself 
the structure of the group under consideration. To go beyond the maximaly-decoupled 
method, we have aimed to construct an SCF theory, i.e., v (external parameter)-dependent 
Hartree-Fock (HF) theory. Toward such an ultimate goal, the v-HF theory has been recon- 
structed on an affine Kac-Moody algebra along the soliton theory, using infinite-dimensional 
fermion. An infinite-dimensional fermion operator is introduced through a Laurent ex- 
pansion of finite-dimensional fermion operators with respect to degrees of freedom of the 
fermions related to a w-dependent potential with a T-periodicity. A bilinear equation for 
the w-HF theory has been transcribed onto the corresponding r-function using the regu- 
lar representation for the group and the Schur-polynomials. The v-BF SCF theory on an 
infinite-dimensional Fock space F^o leads to a dynamics on an infinite-dimensional Grass- 
mannian Groo and may describe more precisely such a dynamics on the group manifold. 
A finite-dimensional Grassmannian is identified with a Groo which is affiliated with the 
group manifold obtained by reducting gl{oo) to sl{N) and su{N). As an illustration we 
will study an infinite-dimensional matrix model extended from the finite-dimensional su{2) 
Lipkin-Meshkov-Glick model which is a famous exactly-solvable model. 
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An original version of this work was first presented by S. Nishiyama at the Sixth International 
Wigner Symposium held in Bogazici University, Istanbul, Turkey, August 16-22, 1999 [1] and 
has been presented by S. Nishiyama at the YITP Workshop on Fundamental Problems and 
Applications of Quantum Field Theory "Topological aspects of quantum field theory" Yukawa 
Institute for Theoretical Physics, Kyoto University, December 14-16, 2006. A preliminary 
version of the recent work has been presented by S. Nishiyama at the KEK String Workshop 
(poster session) 2008 held at High Energy Accelerator Research Organization, KEK, March 
4-6, 2008. 

1 Introduction 

1.1 Historical background on microscopic study of nuclear collective motions 

A standard description of fermion many-body systems starts with the most basic approximation 
that is based on an independent-particle (IP) picture, i.e., self-consistent-field (SCF) approxi- 
mation for the fermions. The Hartree-Fock theory (HFT) is typical one of such an approxima- 
tion for ground states of the systems. Excited states are treated with the random phase ap- 
proximation (RPA). The time dependent Hartree-Fock (TDHF) equation and time dependent 
Hartree-Bogoliubov (TDHB) equation are nonlinear equations owing to their SCF characters 
and may have no unique solution. The HFT and HBT are given by variational method to 
optimize energy expectation value by a Slater determinant (S-det) and an HB wave function, 
respectively [2]. Particle-hole (p-h) operators of the fermions with single-particle states form 
a Lie algebra u{N) [3] and generate a Thouless transformation [3j which induces a representa- 
tion of the corresponding Lie group U{N). The U{N) canonical transformation transforms an 
S-det with M particles to another S-det. Any S-det is obtained by such a transformation of 
a given reference S-det, i.e., Thouless theorem provides an exact wave function of fermion state 
vector which is the generalized coherent state representation (CS rep) of U{N) Lie group [5]. 
Following Yamamura and Kuriyama [6], we give a brief history of methods extracting collec- 
tive motions out of fully parametrized TDHF/TDHB manifolds in SCF. Arvieu and Veneroni, 
and Baranger and independently Marumori have proposed a theory for spherical even nuclei [7] 
called quasi-particle RPA (QRPA) and it has been a standard approximation for the excited 
states of the systems. In nuclei, there exist a short-range correlation and a long-range one [8]. 
The former is induced by a pairing interaction and generates a superconducting state. The 
excited state is classified by a seniority-scheme and described in terms of quasi-particles given 
by the BCS-Bogoliubov theory [9]. The latter is occured by p-h interactions and gives rise to 
collective motions related to a density fluctuation around equilibrium states. The p-h RPA 
(RPA) describes such collective motions like vibrational and rotational motions. It, however, 
stands on a harmonic approximation and should be extended to take some nonlinear effects 
into account. To solve such a problem, the boson expansion HB theory (BEHBT) has been 
developed by Belyaev and Zelevinsky [10], and Marumori, Yamamura and Tokunaga [Tl]. The 
essence of the BEHBT is to express the fermion-pairs in terms of boson operators keeping a pure 
boson-character. The boson representation is constructed to reproduce the Lie algebra of the 
fermion-pairs. The state vector in the fermion Fock space corresponds to the one in the boson 
Fock space by one-to-one mapping. Such a boson representation makes any transition-matrix- 
valued quantity for the boson-state vectors coincide with that for the fermion-state ones. The 
algebra of fermion-pairs and the boson representation have been extensively investigated. The 
fermion-pairs form an algebra so{2N). As for the boson representation, e.g., da Providencia 
and Weneser and Marshalek |12j have proposed boson operators basing on p-h pairs forming 
an algebra su{N). By Fukutome, Yamamura and Nishiyama [13\ I14j . the fermions were found 
to span the algebras so(2iV + l) and so(2A^ + 2) accompanying with u{N + l). The BET ex- 
pressed by Schwinger-type and Dyson-type bosons has been intensively studied by Fukutome 
and Nishiyama [151 HSl IIZl H]- However, the above BET's themselves do not contain any 
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scheme under which cohective degree of freedom can be selected from the whole degrees of 
freedom. 

On the contrary, we have a traditional approach to the microscopic theory of collective motion, 
the TDHF theory (TDHFT) and TDHB theory (TDHBT), e.g., [UllO]. The pioneering idea 
of the TDHFT was suggested by Marumori [7] for the case of small amplitude vibrational 
motions. Using this idea, one can determine the time dependence of any physical quantity, 
e.g., frequency of the small fluctuation around a static HF/HB field. The equation for the 
frequency has the same form as that given by RPA. A quantum energy given by this method 
means an excitation energy of the first excited state. Then, the RPA is a possible quantization 
of the TDHFT/TDHBT in the small amplitude limit. In fact, as was proved by Marshalek and 
Horzwarth [21], the BEHBT is reduced to the TDHBT under the replacement of boson operators 
with classical canonical variables. Using a canonical transformation in a classical mechanics, it is 
expected to obtain a scheme for choosing the collective degree of freedom in the SCF. Historically, 
there was another stream, i.e., an adiabatic perturbation approach. This approach starts from 
an assumption that the speed of collective motion is much slower than that of any other non- 
collective motion. At an early stage of the study of this stream, the adiabatic treatment of the 
TDHFT (ATDHFT) was presented by Thouless and Valatin [22]. Such a theory has a feature 
common to the one of the theory for large-amplitude collective motion. Later the ATDHFT was 
developed mainly by Baranger and Veneroni, Brink, Villars, Goeke and Reinhard, and Mukherjee 
and Pal [23] . The most important point of the ATDHFT by Villars is in introducing a "collective 
path" into a phase space. A collective motion corresponds to a trajectory in the phase space 
which moves along the collective path. Standing on the same spirit, Holtzwarth and Yukawa, 
Rowe and Bassermann [23], gave the TDHFT and Marumori, Maskawa, Sakata and Kuriyama 
so-called "maximally decoupled" method in a canonical form [25]. So, various techniques of 
classical mechanics are useful and then canonical quantization is expected. By solving the 
equation of collective path, one can obtain some corrections to the TDHF result. The TDHFT 
has a possibility to illustrate not only collective modes but also intrinsic modes. However, 
the following three points remain to be solved yet: (i) to determine a microscopic structure of 
collective motion, which may be a superposition of each particle motion, in relation to dynamics 
under consideration (ii) to determine IP motion which should be orthogonal to collective motion 
and (iii) to give a coupling between both the motions. The canonical-formed TDHFT enables 
us to select the collective motion in relation to the dynamics, though it makes no role to take 
IP motion into account, because the TD S-det contains only canonical variables to represent the 
collective motion. Along the same way as the TDHFT, Yamamura and Kuriyama have extended 
the TDHFT to that on a fermion CS constructed on the TD S-det. The CS rep contains not only 
the usual canonical variables but also the Grassmann variables. A classical image of fermions 
can be obtained by regarding the Grassmann variables as canonical ones [26]. The constraints 
governing the variables to remove the overcounted degrees of freedom were decided under the 
physical consideration. Owing to the Dirac's canonical theory for a constrained system, the 
TDHFT was successfully developed for a unified description of collective and IP motions in the 
classical mechanics [27]. 

1.2 Viewpoint of symmetry of evolution equations 

The TDHF/TDHB can be summarized to find optimal coordinate-systems on a group manifold 
basing on Lie algebras of the finite-dimensional fermion-pairs and to describe dynamics on the 
manifold. The boson operators in BET are generators occurring in the coordinate system of tan- 
gent space on the manifold in the fermion Fock space. But the BET's themselves do not contain 
any scheme under which collective degrees of freedom can be selected from the whole degrees of 
freedom. Approaches to collective motions by the TDHFT suggest that the coordinate system 
on which collective motions is describable deeply relates not only to the global symmetry of the 
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finite-dimensional group manifold itself but also to hidden local symmetries, besides the Hamil- 
tonian. Various collective motions may be well understood by taking the local symmetries into 
account. The local symmetries may be closely connected with infinite-dimensional Lie algebras. 
However, there has been little attempts to manifestly understand collective motions in relation 
to the local symmetries. Prom the viewpoint of symmetry of evolution equations, we will study 
the algebro-geometric structures toward a unified understanding of both the collective and IP 
motions. 

The first issue is to investigate fundamental "curvature equations" to extract collective sub- 
manifolds out of the full TDHF/TDHB manifold. We show that the expression in a quasi-particle 
frame (QPF) of the zero-curvature equations described later becomes the nonlinear RPA which 
is the natural extension of the usual RPA. We abbreviate RPA and QRPA to only RPA. We had 
at first started from a question whether soliton equations exist in the TDHF/TDHB manifold 
or do not, in spite of the difference that the solitons are described in terms of infinite degrees of 
freedom and the RPA in terms of finite ones. We had met with the inverse-scattering-transform 
method by AKNS [28] and the differential geometrical approaches on group manifolds [29]. An 
integrable system is explained by the zero-curvature, i.e., integrability condition of connection on 
the corresponding Lie group. Approaches to collective motions had been little from the viewpoint 
of the curvature. If a collective submanifold is a collection of collective paths, an infinitesimal 
condition to transfer a path to another may be nothing but the integrability condition for the 
submanifold with respect to a parameter time t describing a trajectory of an SCF Hamiltonian 
and to other parameters specifying any point on the submanifold. However the trajectory of 
the SCF Hamiltonian is unable to remain on the manifold. Then the curvature may be able 
to work as a criterion of effectiveness of the collective submanifold. From a wide viewpoint of 
symmetry the RPA is extended to any point on the manifold because an equilibrium state which 
we select as a starting point must be equipotent with any other point on the manifold. The 
well-known RPA had been introduced as a linear approximation to treat excited states around 
a ground state (the equilibrium state), which is essentially a harmonic approximation. When 
an amplitude of oscillation becomes larger and then an anharmonicity appears, then we have to 
treat the anharmonicity by taking nonlinear effects in the equation of motion into account. It is 
shown that equations defining the curvature of the collective submanifold becomes fundamental 
equations to treat the anharmonicity. We call them "the formal RPA equation" . It will be useful 
to understand algebro-geometric meanings of large-amplitude collective motions. 

The second issue is to go beyond the perturbative method with respect to the collective 
variables [25]. For this aim, we investigate an interrelation between the SCF method (SCFM) 
extracting collective motions and r-functional method (r-FM) [30] constructing integrable equa- 
tions in solitons. In a soliton theory on a group manifold, transformation groups governing 
solutions for soliton equations become infinite-dimensional Lie groups whose generators of the 
corresponding Lie algebras are expressed as infinite-order differential operators of affine Kac- 
Moody algebras. An infinite-dimensional fermion Fock space -Fqo is realized in terms of a space of 
complex polynomial algebra. The infinite-dimensional fermions are given in terms of the infinite- 
order differential operators and the soliton equation is nothing but the differential equation to 
determine the group orbit of the highest weight vector in the -Fqo [SQ]. The generalized CS 
rep gives a key to elucidate relationship of a HF wave function to a r-function in the soliton 
theory. This has been pointed out first by D'Ariano and Rasetti [31j for an infinite-dimensional 
harmonic electron gas. Standing on their observation, for the SCFM one can give a theoreti- 
cal frame for an integrable sub-dynamics on an abstract F^o- The relation between SCFM in 
finite-dimensional fermions and r-FM in infinite ones, however, has not been investigated be- 
cause dynamical descriptions of fermion systems by them have looked very different manners. 
In the papers [321 [33l [Ml ESI d] , we have first tried to clarify it using SCFM on U{N) group and 
r-FM on that group. To attain this object we will have to solve the following main problems: 
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first, how we embed the finite-dimensional fermion system into a certain infinite one and how 
we rebuilt the TDHFT on it; second, how any algebraic mechanism works behind particle and 
collective motions and how any relation between collective variables and a spectral parameter in 
soliton theory is there; last, how the SCF Hamiltonian selects various subgroup-orbits and how a 
collective submanifold is made from them and further how the submanifold relates to the formal 
RPA. To understand microscopically cooperative phenomena, the concept of collective motion 
is introduced in relation to a TD variation of SC mean-field. IP motion is described in terms of 
particles referring to a stationary mean-field. The variation of a TD mean-field gives rise to cou- 
plings between collective and IP motions and couplings among quantum fluctuations of the TD 
mean-field itself [6], while in r-FM a soliton equation is derived as follows: Consider an infinite- 
dimensional Lie algebra and its representation on a functional space. The group-orbit of the 
highest weight vector becomes an infinite- dimensional Grassmannian Goo- The bilinear equation 
(Pliicker relation) is nothing else than the soliton equation. This means that a solution space of 
the soliton equation corresponds to a group-orbit of the vacuum state. The SCFM does not use 
the Pliicker relation in the context of a bilinear differential equation defining finite- dimensional 
Grassmannian Gm but seems to use implicitly such a relation. In the SCFM a physical concept 
of quasi-particle and vacuum and a coset space is used instead. If we develop a perturbative 
theory for large-amplitude collective motion [25], an infinite-dimensional Lie algebra might been 
necessarily used. The sub-group orbits consisting of several loop-group paths [36] classified by 
the Pliicker relation exist innumerably in Gm so that the SCFM is related to the soliton theo- 
ry in Goo- The Pliicker relation in a coset space ui^m)xU{N-M) E3 becomes analogous with 
the Hirota's bilinear form \38\ I39|. Toward an ultimate goal we aim to reconstruct a theoretical 
frame for a v (external parameter)-dependent SCFM to describe more precisely the dynamics on 
the FoQ. In the abstract fermion Fock space, we find common features in both SCFM and r-FM. 
(i) Each solution space is described as Grassmannian that is group orbit of the corresponding 
vacuum state, (ii) The former may implicitly explain the Pliicker relation not in terms of bi- 
linear differential equations defining but in terms of the physical concept of quasi-particle 
and vacuum and mathematical language of coset space and coset variable. The various BETs 
are built on the Pliicker relation to hold the Grassmannian. The latter asserts that the soliton 
equations are nothing but the bilinear differential equations giving a boson representation of 
the Pliicker relation. The relation, however, has been unsatisfactorily investigated yet within 
the framework of the usual SCFM. We study it and show that both the methods stand on the 
common features, Pliicker relation or bilinear differential equation defining the Grassmannian. 
On the contrary, we observe different points: (i) The former is built on a /inite-dimensional Lie 
algebra but the latter on an in/zniie-dimensional one. (ii) The former has an SCF Hamiltonian 
consisting of a fermion one-body operator, which is derived from a functional derivative of an 
expectation value of a fermion Hamiltonian by a ground-state wave function. The latter intro- 
duces artificially a fermion Hamiltonian of one-body type operator as a boson mapping operator 
from states on fermion Fock space to corresponding ones on r- functional space (r-FS). 

The last issue is, despite a difference due to the dimension of fermions, to aim at obtaining 
a close connection between concept of mean-field potential and gauge of fermions inherent in 
the SCFM and at making a role of a loop group [36] to be clear. Through the observation, we 
construct infinite-dimensional fermion operators from the finite-dimensional ones by Laurent 
expansion with respect to a circle S^. Then with the use of an affine Kac~Moody (KM) algebra 
according to the idea of Dirac's positron theory [30], we rebuilt a TDHFT in Fqo- The TDHFT 
results in a gauge theory of fermions and the collective motion, fluctuation of the mean-field 
potential, appears as the motion of fermion gauges with a common factor. The physical concept 
of the quasi-particle and vacuum in the SCFM on the connects to the "Pliicker relations" due 
to the Dirac theory, in other words, the algebraic mechanism extracting various sub-group orbits 
consisting of loop path out of the full TDHF manifold is just the "Hirota's bilinear form" [39] 
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which is an su{N){£ sl{N)) reduction of gl{N) in the r-FM. As a result, it is shown that an 
infinite-dimensional fermion many-body system is also realizable in a finite-dimensional one 
and that roles of the soliton equation (Pliicker relation) and the TDHF equation are made 
clear. We also understand an SCF dynamics through gauge of interacting infinite-dimensional 
fermions. A bilinear equation for the v-UFT has been transcribed onto the corresponding r- 
function using the regular representation for the group and the Schur polynomials. The f-HF 
SCFM on an infinite-dimensional Fock space F^o leads to a dynamics on an infinite-dimensional 
Grassmannian Groo and may describe more precisely such a dynamics on the group manifold. 
A finite-dimensional Grassmannian is identified with a Groo which is affiliated with the group 
manifold obtained by reducting (7/(00) to sl{N) and su{N). We have given explicit expressions for 
Laurent coefficients of soliton solutions for sl{N) and 'su{N) on the Groo using Chevalley bases 
for sl{N) and su{N) |41j . As an illustration we will attempt to make a f-HFT approach to an 
infinite-dimensional matrix model extended from the finite-dimensional su{2) Lipkin-Meshkov- 
Glick (LMG) model [32]. For this aim, we give an affine KM algebra sl{2,C) (complexification 
of 5m(2)) to which the LMG generators subject, and their t representations and the ax mappings 
for them. We can represent an infinite-dimensional matrix of the LMG Hamiltonian and its HF 
Hamiltonian in terms of the Schur polynomials. Its infinite-dimensional HF operator is also 
given through the mapping 0"^/ for ipiipj of infinite-dimensional fermions ipi and ip* , which is 
expressed by the Schur polynomials Sk{x) and Sk{dx)- Further its r-function for a simple case 
is provided by the Pliicker coordinates and Schur polynomials. 

In Section [21 we propose curvature equations as fundamental equations to extract a collective 
submanifold out of the full TDHB manifold. Basing on these ideas, we construct the curvature 
equations and study the relation between the maximal decoupled method and the curvature 
equations. We further investigate the role of the non-zero curvature arising from the residual 
Hamiltonian. Making use of the expression of the zero-curvature equations in the QPF, we 
find the formal RPA equation. In Section [3l we present a simply unified aspect for the SCFM 
and the r-FM and show a simple idea connecting both the methods. We study the algebraic 
relation between coset coordinate and Pliicker coordinate. Basing on the above idea, we at- 
tempt to rebuilt the TDHFT in r-FS. We introduce u-dependent infinite-dimensional fermion 
operators and a F^o through Laurent expansion with respect to the degrees of freedom of the 
original fermions. The algebraic relation between both the methods is manifestly described. 
We embed a HF u{N) Lie algebra into a gl{oo) by means of infinite-dimensional fermions. 
The u-SCFM in r-FS is developed. The role of the shift operators in the r-FM is studied. 
As an illustration, explicit expressions for Laurent coefficients of soliton solutions for sl{N) 
and 'su{N) are presented. A problem related to a nonlinear Schrodinger equation is also dis- 
cussed. In Section [H we construct a formal RPA equation on F^o and also argue about the 
relation between a loop collective path and a formal RPAEQ. Consequently, it can be proved 
that the usual perturbative method with respect to periodic collective variables in the TDHFT 
is involved in the present method which aims for constructing the TDHFT on the affine KM 
algebra. In Section [5l we introduce infinite-dimensional "particle" and "hole" operators and op- 
erators Kq and K± defined by infinite-dimensional "particle-hole" pair operators. Using these 
operators, we construct an infinite-dimensional Heisenberg subalgebra of the affine KM algebra 
sl{2,C). The LMG Hamiltonian and its HF Hamiltonian are expressed in terms of the Heisen- 
berg basic-elements whose representations are isomorphic to those in the corresponding boson 
space. They are given in terms of infinite numbers of variables and derivatives dx^ through 
the Schur polynomials Sk{x). We give also an infinite-dimensional representation of SU (2A^)oo 
transformation of the particle and hole operators. Finally, in Section [6l we summarize and dis- 
cuss the results and future problems. We give Appendices IAKJI Especially in Appendix [Jj we 
show an explicit expression for Pliicker coordinate for the LMG model and calculate a quantity, 
det(l7v +p^p) {p '■ coset variable) for the LMG model, in terms of the Schur polynomials. 
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2 Integrability conditions and collective submanifolds 
2.1 Introduction 

Let us consider an abstract evolution equation dtu{t) = K[u{t)) for u, which is dependent only 
on a parameter, time t. If there exists a symmetry operation to transfer a solution for u to 
another one, then introducing another parameter s specifying various solutions, we can derive 
another form of evolution equation with respect to s, dsu(t, s) = K(u{t, s)) for which we should 
want to search. The infinitesimal condition for the existence of such a symmetry appears as 
the well-known integrability condition dsK{u{t, s)) = dtK{u{t, s)). The "maximally decoupled" 
method proposed by Marumori et al. [25], invariance principle of Schrodinger equation and 
canonicity condition, can be considered as a theory to apply the above basic idea to certain 
multiple parametrized symmetries. The method is regarded as a mathematical tool to describe 
the symmetries of the collective submanifold in terms of t and collective variables, in which the 
canonicity conditions make the collective variables to be orthogonal coordinate-systems. 

Therefore we adopt a concept of curvature unfamiliar in the conventional TDHBT. The 
reason why we take such a thing is the following: let us consider a description of motions 
of systems on a group manifold. An arbitrary state of the system induced by a transitive 
group action corresponds to any point of the full group parameter space and therefore its time 
evolution is represented by an integral curve in this space. In the whole representation space 
adopted, we assume the existence of 2?7i parameters specifying the proper subspace in which 
the original motion of the system can be approximated well, the existence of the well-defined 
symmetries. Suppose we start from a given point on a space, which consists of t and the 
2m parameters, and end at the same point again along the closed curve. Then we have the 
value of the group parameter different from the one at an initial point on the proper subspace. 
We search for some quantities characterizing the difference of the value. For our aim, we 
introduce a differential geometrical viewpoint. The our basic idea lies in the introduction of 
a sort of Lagrange manner familiar to fluid dynamics to describe collective coordinate systems. 
This manner enables us to take a one-form which is linearly composed of TDHB Hamiltonian 
and infinitesimal generators induced by collective variable differentials of an S0{2N) canonical 
transformation. The integrability conditions of the system read the curvature C = 0. Our 
methods are constructed manifesting themselves the structure of the group under consideration 
to make easy to understand physical characters at any point on the group manifold. 



2.2 Integrability conditions 



We consider many fermion systems with pair correlations. Let Cq, and Ca (a = 1, . . . , A'^) be the 
annihilation-creation operators of the fermion. Owing to the anti-commutation relations among 
them, some sets of fermion operators with simple construction become the basis of a Lie algebra. 



The operators in the fermion so{2N) Lie algebra, E' 



generate a canonical transformation U{g) (the Bogoliubov transformation [9]) which is specified 
by an S0{2N) matrix g: 



[d,d^] = U{g)[c,c^U\g) = [c,c^]g, g = 
U-\g) = U\g), U{g)U{g') = U{gg'), 



a b* 
b a* 



99 



9^9 = l2Af, 



(2.1) 



where (c, c^) = ((cq), (cq)) and (d, d)) = {{di),{d\)) are 2A-dimensional row vectors and a = (a" ) 
and b = (6") (i = 1, . . . ,N) are N x N matrices. I27V is a 2A-dimensional unit matrix. The 
symbols f, * and T mean the hermitian conjugate, the complex conjugation and the transposition, 
respectively. The explicit expressions of the canonical transformations are given by Fukutome 
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for the various types of the fermion Lie algebra [16]. The fermion Lie operators of the quasi- 
particles {E"^ ■, Eij, E"^^) are constructed from the operators d and d) in (|2.1|) by the same way as 
the one to define the set {E'^p, E^p, E'^f^). The set E in the quasi-particle frame is transformed 
into the set E in the particle frame as follows: 



E\ E" 
E.. -E\^ 



9' 



E\ E" 
E.. -E'J 



9- (2.2) 



The E*, = {E^ j) and = (-E*^) etc. are N x N matrices. A TDHB Hamiltonian of the system 
is given by 



' ct " 






c 





-F* -D* 
D F 



T\ (2.3) 



where the HB matrices F = {Fa/3) a-^d D = {D^/s) are related to the quasi-particle vacuum 
expectation values of the Lie operators {E) as 

Fa/3 = Kf3 + [al3\-f6]{{E\) + ^6^s) {F^ = F), 
Daf3 = ^[aj\/36]{Es^) {D'^ = -D), 
[aP\jd] = -[aS\j(3] = Ma(3] = [(3a\6jr. 

The quantities /iq/3 and [a/3|7(5] are the matrix element of the single-particle Hamiltonian and 
the antisymmetrized one of the interaction potential, respectively. Here and hereafter we use 
the dummy index convention to take summation over the repeated index. 

Let |0) be the free-particle vacuum satisfying Cc|0) = 0. The S'0(2A^)(HB) wave function 
\4>{g)) is constructed by a transitive action of the SO{2N) canonical transformation U{g) on |0): 
|<^(^)) = [/-i(^)|0), g G S0{2N). In the conventional TDHBT, the TD wave function \(j){g)) is 
given through that of the TD group parameters a (a*) and b (6*). They characterize the TD 
self-consistent mean HB fields F and D whose dynamical changes induce the collective motions 
of the many fermion systems. As was made in the TDHF case [25] and [27], we introduce a TD 
SO{2N) canonical transformation U{g) = U[g{A{t),A*{t))]. A set of TD complex variables 
(A(t), A*(t)) = (A„(t), A* (t); n = 1, . . . , m) associated with the collective motions specifies the 
group parameters. The number m is assumed to be much smaller than the order of the SO{2N) 
Lie algebra, which means there exist only a few "collective degrees of freedom" . The above is 
the natural extension of the method in TDHF case to the TDHB case. 

However, differing from the above usual manner, we have another way, may be called a Lag- 
range-like manner, to introduce a set of complex variables. This is realized if we regard the 
above-mentioned variables (A(i), A*(i)) as functions of independent variables (A, A*) = (A„, A*) 
and t, where time-independent variables (A, A*) are introduced as local coordinates to specify 
any point of a 2m-dimensional collective submanifold. A collective motion in the 2m-dimensional 
manifold is possibly determined in the usual manner if we could know the explicit forms of A 
and A* in terms of (A, A*) and t. The above manner seems to be very analogous to the Lagrange 
manner in the fluid dynamics. The pair of variables (A, A*) specifies variations of the SCF 
associated with the collective motion described by a pair of collective coordinates a and their 
conjugate vr in the Lagrange-like manner, a = "^(^* + ^) ^ = ~ -^)|25j- Thus, the 

SO{2N) canonical transformation is rewritten as U{g) = U[g{A, A* ,t)] E S0(2N). Notice that 
a functional form g{A{t),A*{t)) changes into another form g{A,A*,t) due to an adoption of the 
Lagrange-like manner. This manner enables us to take a one- form 0, which is linearly composed 
of the infinitesimal generators induced by the time differential and the collective variable ones 
{dt,d\,dA*) of the S0{2N) canonical transformation U[g{A, A* ,t)]. By introducing the one- 
form Q, it is possible to search for the collective path and the collective hamiltonian almost 
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separated from other remaining degrees of freedom of the systems. It may be achieved to study 
the integrabihty conditions of our systems which are expressed as the set of the Lie-algebra- 
valued equations. 

We define the Lie-algebra- valued infinitesimal generators of collective submanifolds as follows: 



H,/h^{idtU-\g))U{g), 
Oi ^ itdA„U-\g))Uig), 



Or. 



{idAM-\g))U{g). 



(2.4) 



Here and hereafter, for simplicity we abbreviate ^(A, A*,t) as g. The explicit form of the 
infinitesimal generators for the TDHF was first given by Yamamura and Kuriyama [27] . In our 
TDHB, the Lie-algebra-valued infinitesimal generators are expressed by the trace form as 



Or,. = 



TV {zdtg • g^) 



-^Tr }^{idA„g . g^] 
iTrl{idA*g.g^] 







E., 




E\ 


E" 


E„ 


-E'J 


E\ 


E" 


E„ 


-E'J 



1 

2 L'-' 



c,c^{idA„g • firl" 



1 

2 L'-' 



c,c^]{idA*g ■ g^ 



(2.5) 



Multiplying the S0{2N) wave function \cl){g)) on the both sides of ([27 
equations on the so(2A^) Lie algebra: 



(2.6) 

we get a set of 



DM9)) = {dt + iH,/hmg))=0 
d 



Da,M9)) = {dA,^+iOimg)) = 0, DaM9)) = (dA^+iOr^mg)) = 0. 



(2.7) 



We regard these equations (|2.7|) as partial differential equations for \(j){g)). In order to discuss 
the conditions under which the differential equations (j2.7p can be solved, the mathematical 
method well known as integrabihty conditions is useful. For this aim, we take a one- form O 
linearly composed of the infinitesimal gener ators ([23]): ^ = -i{Hc/h -dt + Oi- dK + On • c^A* ). 

With the aid of the Q., the integrabihty conditions of the system read C = dO. — Q. f\Q. = 0, 
where d and A denote the exterior differentiation and the exterior product, respectively. From 
the differential geometrical viewpoint, the quantity C means the curvature of a connection. 
Then the integrabihty conditions may be interpreted as the vanishing of the curvature of the 
connection {Dt-, Da,^, E>\^). The detailed structure of the curvature is calculated to be 

C = Ct,A,JK Adt + Ct,A*dK A (it + CA^,,A*dK A dA^' 
+ h^A,,'^ndK A dA„, + iCA*„A*c^A* A dA^,, 



where 



a,A„ = [Dt, DaJ = idtOi - idA„Hj^ + [Ol H,/h], 
Ct,Ai = [Dt, Da{\ = idtOn - idA*„H,/h + [On, Hjh], 
Ca„„K = iDA„„DA^] = idA„,0.n - ^^A^Ol + [0„,Ol,], 
Ca„„a„ ^ [Da„„DaJ = idA^,Oi-idA^Ol + [Olol], 
Ca*,,a* = [Da*,,Da*J = idA*Pn - idA^On' + [0„,0„']. 



{21 



The vanishing of the curvature C means C,^, = 0. 
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Finally with the use of the explicit forms of (j2.5p and (j2.6p . we can get the set of Lie-algebra- 
valued equations as the integrability conditions of partial differential equations 



Cf,A„ = ^[c,C^]Ct,A, 



Ca„,,a„ = ^[c,c"^]Ca„,,a„ 

Ca:„A* = ^[c,ct]CA* ,A* 



Ca,,A* = idAjn - idA*Ol + i^n, Ol], 



(2.9) 



CA„„An = idA„,0i-idAjl + [elel], 

Ca*,,A* = idA* On - idA*On' + [On, On']- 



Here the quantities J-c, On, On are defined through partial differential equations, 
ihdtg = Tc9 and id^^g = O^g, idA*g = Ong- 



(2.10) 



The quantity C,^, may be naturally regarded as the curvature of the connection on the group 
manifold. The reason becomes clear if we take the following procedure quite parallel with the 
above: Starting from (|2.10p . we are led to a set of partial differential equations on the SO{2N) 
Lie group, 



Vtg = {dt + iJ^c/fi)g = 0, 

VA^g = (9a„ + iOi)g = 0, Pa*5 = (9aj; + iOn)g = 0. 



(2.11) 



The curvature C,^, (= [P,,D,]) of the connection (P^jDa^j^^a*) is easily shown to be equiv- 
alent to the quantity C,^, in (j2.9p . The above set of the Lie-algebra- valued equations (|2.9p 
evidently leads us to putting all the curvatures C,^, in (|2.9p equal to zero. On the other hand, 
the TDHB Hamiltonian (j2.3p . being the full Hamiltonian on the full S0{2N) wave function 
space, can be represented in the same form as (12. 4p . HuB/fi = {'idtU~^{g'))U{g'), where g' 
is any point on the SO{2N) group manifold. This Hamiltonian is also transformed into the 
same form as (|2.5p . It is self-evident that the above fact leads us to the well-known TDHBEQ, 
ihdtg' = Tg' . The full TDHB Hamiltonian can be decomposed into two components at the 
reference point g' = g'- 



U-Hg>)=U-Hg) 



Hn + Hj-, 



9' =9 



•Fc~^ -F res ) 



where the second part Hj-cs{J^t-cs) means a residual component out of a well-defined collective 
submanifold for which we should search now. 

For our purpose, let us introduce another curvature Cj and C'^ ^* with the same forms as 
those in ()2.9p . except that the Hamiltonian J^c is replaced by J-\g'=g (= J-'c + J-'rcs)- The quasi- 
particle vacuum expectation values of the Lie-algebra-valued curvatures are easily calculated as 



{Hrcs)c 



-In 
Itv 



g^{J='g - ihdtg ■ g^ 



(2.12) 
(2.13) 



where we have used (Ct,A„) = and (Cf^A*) = 0. The above equations (j2.12p and (|2.13p 
are interpreted that the values of {C't,An)g and (C't^A*)^ represent the gradient of energy of 
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the residual Hamiltonian in the 2m-dimensional manifold. Suppose there exists the well-defined 
collective submanifold. Then it will be not so wrong to deduce the following remarks: the energy 
value of the residual Hamiltonian becomes almost constant on the collective submanifold, i.e., 

6g{H,^s)g = and dA„{H,,,s)g = 0, 5a* (i/rcs)^ = 0, 

where 5g means ^-variation, regarding g as function of (A, A*) and t. It may be achieved if we 
should determine g (collective path) and J-c (collective Hamiltonian) through auxiliary quantity 
{6, 6^) so as to satisfy He + const = -?/hb as far as possible. Putting .Fc = in (j2.9p . we seek 
for g and J-q satisfying 

Ct,A„=0, Ct,A*=0, Ca„„a*=0, Ca„„a„ = 0, Ca;„a*=0. (2.14) 

The set of the equations C,^, = makes an essential role to determine the collective submanifold 
in the TDHBT. The set of the equations (I2.14p and ()2.1ip becomes our fundamental equation 
for describing the collective motions, under the restrictions ()2.2ip . 

If we want to describe the collective motions through the TD complex variables (A(t), A*(t)) 
in the usual manner, we must inevitably know A and A* as functions of (A, A*) and t. For this 
aim, it is necessary to discuss the correspondence of the Lagrange-like manner to the usual one. 

First let us define the Lie-algebra-valued infinitesimal generator of collective submanifolds as 

Ot I {idj^U-\g))U{9), On = {idj,.U~\g))U{9). io G 9), 

whose form is the same as the one in (12. 4p . To guarantee A.„(t) and A*(t) to be canonical, 
according to [25l[27], we set up the following expectation values with use of the S0{2N) (HB) 
wave function \(p{g)): 

{m\i\M9)) = imioum) = ikK, 

imiidA^jm) = imiOnim) = -4^- (2-15) 

The above relation leads us to the weak canonical commutation relation 

m)\[On,dl]\m)=^nn', 

{(t>mololM9)) = 0, {4>mOn,On'M9)) = (n,n' = l,...,m) (2.16) 

the proof of which was shown in [25] and |27j . 

Using (j2.7p . the collective Hamiltonian Hc/h and the infinitesimal generators On and 0„ in 
the Lagrange-like manner are expressed in terms of infinitesimal ones Oi and O n in the usual 
way as follows: 

Hjh = dt'Koi + dtKOn, 

oi = dA^K'Ol + 5a„A^,0„,, On = dA^^An'Ol + ^a* A:,0„,. (2.17) 

Substituting (j2.17p into (j2.8p . it is easy to evaluate the expectation values of the Lie-algebra- 
valued curvatures C,,, by the S0{2N) (HB) wave function \(l)[g{A{t),A*{t))]) (= \(j)[g{A, A* ,t)])). 
A weak integrability condition requiring the expectation values {(j){g)\C,^,\4>{g)) = yields the 
following set of partial differential equations with aid of the quasi-particle vacuum property, 

d\H9)) = 0: 

dA^An'dtK, - dA^KMn' = dM^An'dtK, - dA^K'dtAn' = \i^{n{g)[elF,/h]], (2.18) 
5a*a„»9a„,a;„ -aA^A*„9A„,A„» = \T,{n{g)[en,el]], 
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5a„A„»9a„, A;„ - 9A„A;„aA„, A„" = \T^{n{g)[el ol]}, (2.19) 
5aj,A„»9a;A^„ -aA*A;;„aA;,A„» = iTV{7^(ff)[^„, 

where an S0{2N) (HB) density matrix TZ{g) is defined as 

5^ n\g)=n{g), n\g) = l2N, (2.20) 

in which g becomes function of the complex variables (A, A*) and t. We here have used the 
transformation property (|2.2p . the trace formulae equations (j2.5p and (j2.6p and the differential 
formulae, i.e., 

{^{g)\idf^ Pn\(t){g)) = -\5nn', {(j){g)\idj^* di\^{g)) = ^5nn', 

n n' 

m)\id^ oMg)) = 0, imiidA.pnm)) = o, 

which owe to the canonicity condition (j2.15p and weak canonical commutation relation (12.16p . 

Through the above procedure, as a final goal, we get the correspondence of the Lagrange- 
like manner to the usual one. We have no unknown quantities in the r.h.s. of equations (j2.18p 
and (j2.19p . if we could completely solve our fundamental equations to describe the collective 
motion. Then we come up to be able to know in principle the explicit forms of (A, A*) in terms 
of (A, A*) and t by solving the partial differential equations (I2.18P and (|2.19p . However we 
should take enough notice of roles different from each other made by equations (|2.18p and (j2.19p , 
respectively, to construct the solutions. Especially, it turns out that the l.h.s. in (j2.19p has a close 
connection with Lagrange bracket. From the outset we have set up the canonicity condition to 
guarantee the complex variables (A, A*) in the usual manner to be canonical. Thus the variables 
(A, A*) are interpreted as functions giving a canonical transformation from (A, A*) to another 
complex variables (A, A*) in the Lagrange-like manner. From this interpretation, we see that 
the canonical invariance requirements impose the following restrictions on the r.h.s. of p.l9p : 

-lTr{n{g)[Bn,9l]} = dnn', lTr{n{g)[0lel]} = O, iTV{7^(5)[^n, ^n']} = 0. (2.21) 

Using (j2.19p and ()2.2ip . we get Lagrange brackets for canonical transformation of (A, A*) 
to (A,A^). 

2.3 Validity of maximally-decoupled theory 

First we transform the set of the fundamental equations in the particle frame into the one in 
the quasi-particle frame. The SO{2N) (TDHB) Hamiltonian of the system is expressed as 

H = To, (2.22) 

the relation of which to the original TDHB Hamiltonian T is given by To = g'^Tg, g G SO{2N). 
The infinitesimal generators of collective submanifolds and their integrability conditions ex- 
pressed as the Lie-algebra- valued equations are also rewritten into the ones in the quasi-particle 
frame as follows: 

H, = ^[d,cl)]To-c 



-In 
Itv 



H^^ = \[d,S]To 



d 



To 



-F* 

o 



-D* 

^ o 



d 



On = \[d,S]eo- 



d 
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Ct,A„ = ^[d,d^]Co^ 



t,An 



d 



0, 



Ct,Kl = \[d,S]Co- 



0, 



d 

Ca„,,A„ = ^['i,C?"^]Co-A„,,A„ ^ =0, 

Co-A„,,A„ = idK^,B\_^ - id^jl^^, 
Ca„,,A* = 5['i,C^'^]Co-A„,,A* ^ 
Co-A„,,A* = idA^,do~n - idA*el_„, 

Ca;,,a* = |[(i,d"^]Co-A;„A* 

Co- A*,, A* = idA^Oo-n - idAlOo-n 



n O—W 

d 



b'^ 1 

'o~rii o—n'\i 



(2.23) 



The quantities J^o-c, ^l-n (— d^^ng) and 6o-n (= g^^^ng) are defined through partial differential 
equations on the SO{2N) Lie group manifold, 



ihdtg'' = To^cg 



and 



idAtg^ = Oo-ng 



(2.24) 



In the above set of ()2.23|) . all the curvatures Co-»,» should be made equal to zero. 

The full TDHB Hamiltonian is decomposed into the collective one and the residual one as 



-f^HB — Hc + Hj-cs, To — To-c + To- 



(2.25) 



at the reference point g on the S0{2N) group manifold. Following the preceding section, let us 
introduce other curvatures and C[ with the same forms as those in (I2.23p except that 

To-c is replaced by Tq. Then the corresponding curvatures C'^_t and C^_jy\^* are also divided 
into two terms, 



/ic I /ires 1^1 /-ic I /ire 

•-O-t.An ~ ^0-t,An ^ ^0-t,An^ ^0~t,A* " ^0-t,A* ^O- 



Here the collective curvatures C 



o-t,A„ 



and C 



o~t,A* 



arising from To-c are defined as the same 
forms as the ones in (12.23p . The residual curvatures C^^j and C^^^ ^* arising from To-rcs are 
defined as 

^o-tM = -idAr.To^rcs/h- [el^n, Torres 

Co-t,A* = -idA*^To-rcs/h - [Oo-n, ^o-rcs/^]- (2.26) 

Using (|2.24p and (I2.25P , the Lie-algebra- valued forms of the curvatures are calculated as 



t,An 



-idA„Hrcs/h, C\ 



■res 
t,A* 



-idAtHxcsl^- 



Supposing there exist the well-defined collective submanifolds satisfying ()2.24p . we should de- 
mand that the following curvatures are made equal to zero: 



Co-t,A„ — 0' '^o-t,A* — 0) 

Co-A„,,A„ = 0, Co_A„,,A* = 0, Co_A* ,A* = 0, 



(2.27) 
(2.28) 
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the first equation (j2.27p of which lead us to the Lie- algebra- valued relations, 



C' 



t,An 



/^res 



-idA*^Hrcs/fi- 



(2.29) 



Then the curvature C'^j^^^ and C'^j^^* can be regarded as the gradients of quantum-mechanical 
potentials due to the existence of the residual Hamiltonian Hj-^s on the collective submanifolds. 
The potentials become almost flat on the collective submanifolds, i.e., i^HB = -ffc + const, if the 
proper subspace determined is an almost invariant subspace of the full TDHB Hamiltonian. This 
collective subspace is an almost degenerate eigenspace of the residual Hamiltonian. Therefore 
it is naturally deduced that, provided there exists the well-defined collective subspace, the 
residual curvatures at a point on the subspace are extremely small. Thus, the way of extracting 
the collective submanifolds out of the full TDHB manifold is made possible by the minimization 
of the residual curvature, for which a deep insight into (j2.29p becomes necessary. 

Finally, the restrictions to assure the Lagrange bracket for the usual collective variables and 
Lagrange-like ones are transformed into the following forms represented in the QPF: 



-In 
Itv 

-In 
In 



[0^ 



o-n, ^o-n'i 



■n' o—n' 



0, 



Tr 



-Itv 
Itv 



^o— n'] 



We discuss here how the Lagrange-like manner picture is transformed into the usual one. 
First let us regard any point on the collective submanifold as a set of initial points (initial value) 
in the usual manner. Suppose we observe the time evolution of the system with various initial 
values. Then we have the following relations which make a connection between the Lagrange-like 
manner and the usual one 



:Fo-c/h = dtAJl-n + dtK^o-n, 



(2.30) 
(2.31) 



in which the transformation functions are set up by the initial conditions, 



A„(t)|t=o = A„(A, A*,t)|t=o = A„, 



A*(^)|,=o = A^(A,A^^)|i=o = A*, 

5nn', 9A„A*/|i=o = 0, 9A*A„/|f=o 



0, 



in order to guarantee both pictures to coincide at time t = Q. On the other hand, our collective 
Hamiltonian J-q-c can also be expressed in the form 



To^Jh = 7;„(A, K\t)el_^ + <(A, k\t)eo-n, 



(2.32) 



where the expansion coefficients Vn and are interpreted as velocity fields in the Lagrange-like 
manner. Substituting (j2.3ip into (j2.32p and comparing with (|2.30p . we can get the relations 



A„ = dtkn = Vn'dK An + <,9a*, A„, = A* = Vn'dA + v*,dA*^Al, 

from which the initial conditions of the velocity fields are given as 

A„(t)|t=o = dtAn\t=o = ^;„(A,A^^)|t=o, kiit)\t=o = 9jA*|t=o = <(A, A^ t)|t=o- 



Then we obtain the correspondence of the time derivatives of the collective co-ordinates in the 
usual manner to the velocity fields in the Lagrange-like one. 
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Finally we impose the canonicity conditions in the usual manner, 

{m\oi\m) = ihK, {m\On\m) = -^Vk, (2.33) 

which leads us to the weak canonical commutation relation with the aid of (|2.28|) and (j2.2ip . 

The TDHBT for maximally- decoupled collective motions can be formulated parallel with TD- 
HFT [25] . The basic concept of the theory lies in an introduction of the invariance principle of 
the Schrddinger equation, and the TDHBEQ is solved under the canonicity condition and the 
vanishing of non-collective dangerous terms. However, as we have no justification on the valid- 
ity of the maximally- decoupled method, we must give a criterion how it extracts the collective 
submanifold effectively out of the full TDHB manifold. We are now in a position to derive some 
quantities by which the criterion is established. For this aim, we express the collective Hamilto- 
nian J-q-c and the residual one J-o-res in the same form as the one of the TDHB Hamiltonian To 
given in ()2.22p . We also represent quantities 0\_^, ^o-t A„ ^^'^ ^o-t A* which consist of N x N 
block matrices as follows: 



ni - 


Co ^0 


- n 




T 

'^0 = 




(2.34) 


/ires 




L '-^ 


/-•res 

/^resT 
-4 J 


n 


pres 


presT /ires 




^res 
'~'o~t,A*„ 




L 


nrcs 

J 


pres 

Lg* - - 

n 


pres] 


rircs /-ires* 


C-^ = C--. (2.35) 



Substitution of the explicit form of J^o-rcs and equations p. 340 and p.35p into ()2.26p yields 

= -idA^^Do-res/h + ClnDo^rcs/h + Do^res / hCo,n + ^o,nF*_,^J H + Fo-rcs/Ho,n, (2.36) 
^ip,n ~ ^^A„ -Co-res/^ ~^ Co,nDo-T:cs/^ + -^o-rcs/^Co,n ~ 'fo,nFo-rcs / ^ ~ Fo-resl ^'^o,n- 

The quantity ^o_„ can also be expressed in the same form as the one defined in (|2.34p . Substi- 
tuting this expression into ()2.30p and (j2.33p . we obtain the relations 

-Fo-res//i = Fo/h + dtKiln + dtKC,n. 

Do-rcs/h = Do/h + dtAni>ln + dtK€,n, (2-37) 

together with their complex conjugate and 

Ttio,n = iK^ Tr|t„ = _iA„, (2.38) 

where we have used ()2.25p and the explicit forms of the Hamiltonian. 

As was mentioned, the way of extracting collective submanifolds out of the full TDHB mani- 
fold is made possible by minimization of the residual curvature. This is achieved if we require 
at least expectation values of the residual curvatures to be minimized as much as possible, i.e., 

iHgWtJMg)) = iTVq- - 0, {H9)\Ci:i.\m) = ^TrC|f„ - 0. (2.39) 

We here adopt a condition similar to one of the stationary HB method as was done in the 
TDHF [27]: The so-called dan gerous terms in the residual Hamiltonian To— res 

are made to 

vanish. 



Do— res — 0) 



D* = 0. 

o— res 



(2.40) 
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With aid of equations ([236]) . ([2371) and (p^Sl) . equations (ITIO]) and ([239]) are rewritten as 

Do/h = -dtK^l^ - dtKifln, Dl/h = -dtKipo,n - OtK^Pln, (2.41) 
idA,jT Forres /n = idA^TiFo/h - 5a„ {dtkn'K' - dtK'K') = 0, 

idA^Tr Forres /h = idA*TiFo/h - 5a* {dtK'K, - dtK'^^n') = 0. (2.42) 



First, we wih discuss how equation (j2.4ip leads us to the equation of path for the cohective 
motion. Notice that the quantities 0\_^ and 6o-n are subjected to satisfy the same type of partial 
differential equation as that of (j2.24p . Remember the explicit representation of an SO{2N) 
matrix g given in the previous section. Then we have partial differential equations 

^o,n = -i{dA„b^a + d^yh), (^o,n = -i(5A„at6'^ + dj,J^d*). (2.43) 

together with its complex conjugate. Putting the relation = g'^ ^9 and (j2.43p into (j2.4ip . we 
get 

{ {Da + Fb)/h- {iKd'^^ b + i A^^^. h) } 

+ b^{- {F*d + D*b)/h - {iknd^^d + iA^a^.a)} = 0, (2.44) 

Let H be an exact Hamiltonian of the system with certain two-body interaction and let us 
denote the expectation value of H by \4>ig)) as {H)g. It can be easily proved that the relations 

da*{H)g = -l{F''a + D%), db*{H)g = ^{Fb + Da), (2.45) 

and their complex conjugate relations do hold, through which the well-known TDHBEQ is 
converted into a matrix form as 



{H)g/^/n. (2.46) 



The quantity {H) means now in turn an expectation value of being a function of 

a/\/2, b/^/2 and their complex conjugate. With the aid of a relation similar to p.45p . equa- 
tion (j2.44p is reduced to 



+ [d,.,V2W„^Jh - + ^kdKt^) I = 0, (2.47) 

As one way of satisfying (j2.47p . we may adopt the following type of partial differential equations: 

5aVv^^^)s/V2/^- (^^n^A„«/v^ + ^A;9A.a/V2) = 0, 

^^lV2^H)-,^Jh - {iLd^J/V2 + iA^^A J/^/2) = 0. (2.48) 

Here we notice the invariance principle of the Schrddinger equation and the canonicity condition 
which leads us necessarily to the equation of collective motion expressed in the canonical forms 



iAt 



-dK (^)§/V2/^' = {H)^^/^JK (2.49) 



which can be easily derived with the use of (|2.46p and (|2.38p . Instead of solving approximately 
our nonlinear time evolution equation (j2.27p . we adopt the above canonical equation. Then, as 
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is clear from the structure of (j2.48p . it is self-evident that equation (|2.48p becomes the equation 
of path for the collective motion under substitution of (|2.49|) . In this sense, equation (|2.48|) is the 
natural extension of the equation of path in the TDHF case |25^I27| to the one in the TDHB case. 
The set of (|2.48p and (|2.49p is expected to determine the behaviour of the maximally decoupled 
collective motions in the TDHB case. However, it means nothing else than the rewriting of 
the TDHBEQ with the use of canonicity condition, if we are able to assume only the existence 
of invariant subspace in the full TDHB solution space. The above interpretation is due to the 
natural consequence of the maximally decoupled theory because there exists, as a matter of case, 
the invariant subspace, if the invariance principle of the Schrddinger equation does hold true. 
The maximally decoupled equation can be solved with the additional RPA boundary condition, 
though its solution is, strictly speaking, different from the true motion of the system on the full 
SO{2N) group manifold. But how can we convince that the solution describes the well-defined 
maximally decoupled collective motions from the other remaining degrees of freedom of motion? 
Therefore, in order to answer such a question, we must establish a criterion how we extract the 
collective submanifolds effectively out of the full TDHB manifold. 

Up to the present stage, equation (j2.42p remains unused yet and makes no role for approa- 
ching to our aim. Finally with the aid of (I2.42p . we will derive some quantity by which the 
range of the validity of the maximally decoupled theory can be evaluated. As was mentioned 
previously, we demanded that the expectation values of the residual curvatures are minimized as 
far as possible and adopted the canonical equation in place of our fundamental equation ()2.27p . 
Then, by combining both the above propositions, it may be expected that we can reach our final 
goal of the present task. Further substitution of the equation of motion ()2.49p (rewrite (H)^^^ 
in the original form {H). again) into (j2.42p yields 



5A„A„'Tr 



(di n-kind\, . n-An"dl . n)^ 



+ 5A„A^.Tr 



n n n it ' 1 1, 



n" ^ n' ft 



(2.50) 



Here we have used the transformation property of the differential 9a„ = 9a„A„/9^ ^ -|-9a„ A*,9y^* 
and the differential formulae for the expectation values of the Hamiltonians H and i^HB 



2.4 Nonlinear RPA theory arising from zero-curvature equation 

Our fundamental equation may work well in the large scale beyond the RPA as the small- 
amplitude limit. A linearly approximate solution of the TDHBEQ becomes the RPAEQ. Suppose 
we solve the fundamental equation by expanding it in the form of a power series of the collective 
variables A and A*(n = 1, . . . ,m; m <C N{2N — l)/2) defined in the Lagrange-like manner. 
Then we must show that the fundamental equation has necessarily the RPA solution at the 
lowest power of the collective variables which approach in the small amplitude limit. A paired 
mode amplitude g{A., A*,t) is separated into stationary and fluctuating components as g = g^°^g. 
This means that the SO{2N) matrix g is decomposed into a product of stationary matrix g^"^ 
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and g{A,A*,t) (~ g). The stationary g^°^ satisfies the usual static S'0(2iV)(HB) eigenvalue 
equation. 

Using the above decomposition of g, an original 50(2A^)(HB) density matrix TZ{A, A*, t) and 
a HB matrix !F{A, A* ,t) are decomposed as 7^ = g^°^'R,g^°^^ and !F = g^°^J-g^''^\ respectively. 
The fluctuating TZ and the HB matrix in fluctuating QPF are given in the following forms: 



2R{g) - In 



-2K*ig) 
-2R*{g) + 1^ 



(2.51) 
(2.52) 



in which all the quantities are redefined in [T3]. Quasi-particle energies e^"^ include a chemical 
potential. 

Introducing fluctuating auxiliary quantities On = g^°^^Ong^°^ and On = g^"^^Ong^°\ then under 
the decomposition g = g^°^g, the zero-curvature equation C, , = in (j2.9p is transformed to 



idtOi - idK^rjh + [olrjh] = o, idtOn - idK*^T,in + \0n. T,in\ = o, 

idKjn-idK*Jl + [0n,0U=^^ 

i9A„, Oi - idAjl + i^l = 0' id^.^On - idA^Jn' + On'] = 0, 



(2.53) 



(2.54) 



lTv{n{g)[On, Ol]} = 5nn', iTr{7^(^) [^1 , Ol]} = 0, lTr{n{g)[On, 4']} = 0, (2.55) 



where the quantities J-c, On and On satisfy partial differential equations, 
ihdtg = Teg, idx^g = Olg and idA*J = Ong. 



(2.56) 



Putting Tc = J' ()2.52p in (|2.53p . we are able to look for a collective path (^) and a collec- 
tive Hamiltonian (J^c) under the minimization of the residual curvature arising from a residual 
Hamiltonian [Txe^- Next, for convenience of further discussion, we introduce modified fluctua- 
ting auxiliary quantities 0\_^ = g^Ong and 6'o-n = 9^ Ong- Then we can rewrite our fundamental 
equations (12.53p . (12.54p and (j2.55p in terms of the above quantities as follows: 



i'dtOl_n 



i~g'^ [dA^TM~g = ^, 



idtOo- 



ig^ {dA*:Fjh)g = o 



(2.57) 



idA,Jo-n - idA*jl_n, " [^o-n, O^n'] = 0, ^QaJI 
idA* Oo~n — idA^Oo-n' — [Oo-n,0o-n'] = 0, 



idAjl_n' 



(2.58) 



iTr 

4 



-In 
Itv 

-In 
In 



'o-ni "o-n'i 



0, 



iTr 

4 



-In 
In 



[Oo^n , Oo^n'] } = 0. (2.59) 



In the derivation of equations (I2.57P and ()2.58p . we have used ()2.56p . The equation (I2.59P is 
easily obtained with the aid of another expression for the fluctuating density matrix TZ (j2.5ip . 



na) = 9 



-In 
Itv 



it 



(2.60) 



In order to investigate the set of the matrix- valued nonlinear time evolution equation (I2.57P 
arising from the zero curvature equation, we give here the dA„ and Oa* differential forms of the 
TDHB density matrix and collective hamiltonian. First, using p.60p . we have 



5An^(5) = dAj 



-In 
In 



9^ + 9 



-In 
Itv 
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-In 
In 



a^ +9 



-In 
In 



dAjHaa^), 



(2.61) 



where we have used the relation g'^g = gg'^ = 1. Using (j2.56p . the above equation is written as 

L Iat J [ u Iat J 

= -^{9m~9 \ -\ . N + ^4 /I 

[ U Iat J [ U iiv J 

Next, by using 6'o_„ = g^Ong and 6o-n = g'^^ng, equation (j2.6ip is transformed into 



dA„T^ig) = -ia ^ 5^ 

L ^ J _ 

On the other hand, from (|2.51|) . we easily obtain 

dAr^ng) 



{2.62) 



2dA^R{g) -2dA„K%g) 
2dAj<{~g) -2dA^R^{~g) 



(2.63) 



Let us substitute explicit representations for g 



into the r.h.s. of (j2.62|) and combine it with (|2.63|) . Then, we obtain the final c?a„ differential 



a h* 
b a* 



and 9 



'4'o-n ~Co— n 



form of the TDHB {SO{2N)) density matrix as follows: 

dAr,R{g) = i{b*ipo~na' - aipo^nb^), dA„R*{g) = -i{a*tpo-J>^ - bipo-^nO^), 
dA^K{g) = i[d*ipo-na^ - bifo-J)^), dA„K*{g) = -i(b*'ipo-nb'^ - afo-nO^)- 

The (9a* differentiation of the S0{2N) density matrix is also made analogously to the above. 

As shown in [H], the fluctuating components of the HB matrix T ()2.52p are linear functionals 
of R{g) and K[g). We can easily calculate the Ba^ differential as follows: 

dAj = {■D)dA„K{~g) + (D)dA^K*{~g) + {d)dA^R{~g) = i{B}i;o-n + i(D}^o-n, 
dA„r = {F)dA„K{g) + iT)dA„k*{g) + indA^m = i{F*}^o-n + i{T}ipo-n. 

Here the matrices (D) etc. are given in [Hj. Similarly, new matrices (D} etc. are defined by 

{B} = \\{ij\D\kl}\\, {B} = \\iij\D\kl}\\ and 

(F*} = \mF*\kl}\l {T} = \\mT\kl}\\, (2.64) 
iij\D\kl} = iij\D\k'iral,kat'i - {ij\D\k'l')bl,,bt,i + iij\d\k'l')bl,,dt,i, 
-{ij\D\kl} = {ii\D\k'l')bk'kh'i - {ij\D\k'l')dk,kdi,i + {ij\d\k'l')dk'kbi'l, 
{ij\F^\kl} = iij\F*\k'l')dl,,dt,i - iij\T\k'l')bl,,^i + {ij\nk'l'%,d^,i, 
-{ij\T\kl} = {ij\F*\k'l')h,kbi'i - {ij\T\k'l')dk'kai'i + {ij\nk'l')dk'kbi'i. 



The above summation is made with indices k' and /' (1 ^ N). Putting JF^, = JF, we have 



dA„^c = i 



(D}Vo-n + (D}Vo-n, ^{F*}iJo~n +^ {F^}^o-n 



(2.65) 



where (F*} etc. stand for the matrices in which the indices i and j in (|2.64p are exchanged. 
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We here derive a new equation formally analogous to the SO{2N) RPA equation. To 
achieve this, we first further decompose the fluctuating pair mode amplitude g into a prod- 
uct of a fluctuating S0{2N) matrix and a 2A^-dimensional diagonal matrix with an exponential 
time dependence as follows: 



9 99i£, -e), 

£ — (5jj£'j) , 



exp [iet /h], 

0, exp[—iet/h] 



(2.66) 



where we redenote a new fluctuating pair mode as g and ej is the A and A* dependent quasi- 
particle energy including the chemical potential. Next, using (j2.56p . = g'^Ong and 6o-n = 
g^On.g, the modified fluctuating auxiliary quantities 6o_n can be written as 



9^ 

'o—n 



g^e,-e)9l_J{e,-e) + 



-dA^et/h, 
0, dA„£t/h 



where we again redenote the new fluctuating auxiliary quantities as 0o_„. Accompanying the 
above change, dtO\_^ are modified to the following forms by using the explicit expression for 0\_^: 



dtio-n - dA^e/h - i[£/h,^o-n], dt^fo-n " ii^/K iPo-n] + 



(2.67) 
g(e,-e). 



In the above, hereafter we adopt (A, A*)-independent e^°^ given in (j2.52p as the quasi-particle 
energy e. If we substitute equations (j2.65p . (|2.66p and ()2.67p into the set of the matrix- valued 
nonlinear time evolution equation, i.e., the equation of (j2.57p . we finally obtain the following 
set of matrix-valued equations: 



o—n J 



. +{-D}i;o-n + {'D}^o~n +T{F*}V^o_„ + T{F}v9o_ 



X g{£ 



(0) _£{0) 



(2.68) 



with the modified new matrices {D} etc. defined through 

{D} = \\{ij\D\kl}\\, {D} = \\{ij\D\kl}\\ and 
{F*} = ||fe1i^1fc/}||, {r} = ||{^j|r|A:/}||, 

whose matrix elements are given by 

{ij\D\kl} = ai>idj'j{i'j'\D\kl} - bi'ibj'j{i'j'\D\kl}* 

- Wiafj{i'j'\F*\kl} + di>ibj'j{i'j'\T\kl}\ 
-{ij\D\kl} = bi,ibj,j{i'j'\D\klY - di,idj,j{i'j'\D\kl} 

- di'ibjij{i'j'\F*\klY + biiidjij{ij'\F*\kl}, 
-{ij\F*\kl} = b*,idj,j{'ij'\D\kl} - d*,,fij,j{ij'\D\klY 

- dUdf,{i'f\F''\ki} + bt,,b,,,ii:j'\T\kiY 

{ij\T\kl} = d^,,bj,,{i'j'\D\klY - ^,ray^{i'j'\D\kl} 
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- b*,,byjii'f\F*\kir + dt,fa,,j{i'j'\F*\kl}, 

in which summation is made over indices i' and j' running from 1 to N. In (j2.68p by making 
block off-diagonal matrices vanish, we get a TD equation with respect to ipo-n and ^Po-n which 
is formally analogous to that of the SO{2N) RPA, though our TD amplitude and matrices {D} 
etc. have (A, A*, t)-dependence. 



2.5 Summary and discussions 

We have studied integrability conditions of the TDHBEQ to determine collective submanifolds 
from the group theoretical viewpoint. As we have seen above, the basic idea lies in the introduc- 
tion of the Lagrange- like manner to describe the collective coordinates. It should be noted that 
the variables are nothing but the parameters to describe the symmetry of TDHBEQ. By intro- 
ducing the one- form, we gave the integrability conditions, the vanishing of the curvatures of the 
connection, expressed as the Lie-algebra-valued equations. The full TDHB Hamiltonian /^hb is 
decomposed into the collective Hamiltonian He and the residual one -ffrcs- To search for the well- 
defined collective submanifold, we have demanded that the expectation value of the curvature is 
minimized so as to satisfy iifres — const or He + const = Hub as far as possible. Further we have 
imposed the restriction to assure the Lagrange bracket for the usual variables and Lagrange-like 
ones. Our fundamental equation together with the restricted condition describes the collective 
motion of the system. 

We have proposed the minimization of the residual curvature arising from the residual part of 
the full TDHB Hamiltonian to determine the collective submanifold. With our theory it is also 
possible to investigate the range of the validity of the maximally decoupled theory of the TDHBT 
with use of the condition to satisfy p.50p . This condition makes an essential role to give the 
criterion how we extract well the collective submanifold out of the full TDHB manifold. The 
reason why the condition occurs in our theory which did not appear in the maximally decoupled 
theory lies in the consideration of the d'^d-type in the residual Hamiltonian to calculate the 
residual curvature and in the adoption of the canonical equation. Since the maximally decoupled 
theory has no consideration of such type from the outset, the condition is trivially fulfilled. This 
is the essential difference between the maximally decoupled theory and ours. 

We have investigated the nonlinear time-evolution equation arising from zero-curvature equa- 
tion on TDHB {SO{2N) Lie group) manifold. It is self-evident that the new equation has an 
SO{2N) RPA solution as a small-amplitude limit. The new equation depends on the collective 
variables (A, A*) defined in a Lagrange-like manner. It works well in the large scale beyond the 
SO{2N) RPA under appropriate boundary and initial conditions. The integrability condition 
is just the infinitesimal condition to transfer a solution to another solution for the evolution 
equation under consideration. The usual treatment of the RPA for small amplitude around 
ground state is nothing but a method of determining an infinitesimal transformation of sym- 
metry under the assumption that fluctuating fields are composed of only normal-modes. We 
conclude that the set of equations defining the symmetry of the SCF equation and the weak 
boson commutation relations on the QPF becomes the nonlinear RPA theory. 

Finally, following Rajeev [l3], we also show the existence of the homogeneous symplectic 
2-form uj. From (12:20]) . using the coset variable q (= fea"!) = -q^, the SO{2N) (HB) 

density matrix TZ{g) is expressed as 



na) = 9 



-In 
1^ 



2R{g) - In -2K*{g) 
2K{g) -2R*{g) + 1^ 



R{9)=q'q{lN + qU) \ K{g)=-q{lN + qU) (2.69) 
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± 

Introducine a new 

x2iV matrix as Z(5r) = (lAT + gtg) 2 [1^, gt] ^ith 2:(5)tZ(5) = I27V, 
we have a very simple expression for TZ{g) as 

n{g) = l2N -2Z{g)Z{gy = -2 \ ^ . ^ ^ , 

Tr7^(5) =0, 

which has quite the same form as the one given by Rajeev [l3j. The two- form lo is given as 



CO 



YTt{ {dn{g)f } = -|Tr{ {dn{g)f n{gf}, 



duj = —duj = (closed form). 
If we introduce hermitian matrices U 



u 
v) 



and V 



V 
t;t 



then we have 



L0{U,V) 



which is a symplectic form and makes it possible to discuss geometric quantization on a finite/in- 
finite-dimensional Grassmannian [44l 1451. 



3 SCF method and r-functional method on group manifolds 
3.1 Introduction 

Despite the difference due to the dimension of fermions mentioned in Section [H we ask the follo- 
wing: How is a collective submanifold, truncated through the SCF equation, related to a subgroup 
orbit in the infinite-dimensional Grassmannian by the r-FM? To get a microscopic understan- 
ding of cooperative phenomena, the concept of collective motion is introduced in relation to TD 
variation of a SCF. Independent-particle (IP) motion is described in terms of particles referring 
to a stationary MF. The TD variation of the TD SCF is attributed to couplings between the 
collective and the IP motions and couplings among quantal fluctuations of the TD SCF [6]. 
There is a one-to-one correspondence between MF potentials and vacuum states of the system. 
Decoupling of collective motion out of full-parameterized TDHF dynamics corresponds to trun- 
cation of the integrable sub- dynamics from a full-parameterized TDHF manifold. The collective 
submanifold is a collection of collective paths developed by the SCF equation. The collectivity 
of each path reflects the geometrical attribute of the Grassmannian, which is independent of the 
characteristic of the SCF Hamiltonian. Then the collective submanifold should be understood 
in relation to the collectivity of various subgroup orbits in the Grassmannian. The collectivity 
arises through interference among interacting fermions and links with the concept of the MF 
potential. The perturbative method has been considered to be useful to describe the periodic 
collective motion with large amplitude |25^ [6] . If we do not break the group structure of the 
Grassmannian in the perturbative method, the loop group may work under that treatment. 

Thus we notice the following point in both methods: Various subgroup orbits consisting of loop 
path may infinitely exist in the full-parameterized TDHF manifold. They must satisfy an infinite 
set of Pliicker relations to hold the Grassmannian. As a result, the finite-dimensional Grass- 
mannian on the circle is identified with an infinite-dimensional one. Namely the r-FM works 
as an algebraic tool to classify the subgroup orbits. The SCF Hamiltonian is able to exist in the 
infinite-dimensional Grassmannian. Then the SCFT can be rebuilt on the infinite-dimensional 
fermion Fock space and also on the r-functional space. The infinite-dimensional fermions are 
introduced through Laurent expansion of the finite-dimensional fermions with respect to the 
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degrees of freedom of the fermions related to the MF potential. Inversely, the collectivity of 
the MF potential is attributed to gauges of interacting infinite-dimensional fermions and inter- 
ference among fermions is elucidated via the Laurent parameter. These are described with the 
use of affine KM algebra according to the Dirac theory [lO]. Algebro-geometric structure of 
m/iniie-dimensional fermion many-body systems is realized in the /inite-dimensional ones. 



3.2 Bilinear differential equation in SCF method 



Owing to the anti-commutation relations {cQ,,ct} = 6a(3, {ca,cp} = {ca,ct} = 0, fermion pair 



operators 6^/3 = c^C/j satisfy a Lie commutation relation [cap, e^s] = Sp^CaS — Sase^f^ and span 

the u{N) Lie algebra. A canonical transformation U{g) = e^°'i^'^°''^i^ (7^ = —7, g = 3 U{N)) 
generates a transformation such that 

U{g)c^aU-^{g) = c'^g/sa, U{g)caU-^{g) = cpg^^, 

U'\g) = U{g~^) = U{g^), U{gg') = U{g)U{g'), g^g = gg^ = lN. (3.1) 
Let |0) be a free vacuum and \4>ai) be an M particle S-det 
c„|0)=0, a = l,...,N, = cir--c\\0), 

U{gMM) = {c^9)m ■ ■ ■ (cMilO) = \g), U{gm = |0), (3.2) 

where the is an iV-dimensional row vector — (^l? * * * ? '-'jv)* Equation (j3.2p shows that the M 
particle S-det is an exterior product of M single-particle states and that U{g) transforms |</>m) 
to another S-det (Thouless transformation) [1] under (j3.ip . Such states are called "simple" 
states. The set of all the "simple" states of unit modulus together with the equivalence relation, 
identifying distinct states only in phases with the same state, constitutes a manifold known 
as Grassmannian Gtm- The Gtm is an orbit of the group given through (13. 2p . Any simple 
state \4>m) G Gtm defines a decomposition of single-particle Hilbert spaces into sub-Hilbert 
spaces of occupied and unoccupied states [46]. Thus, the Gvm corresponds to a coset space 



Grjv/ ~ (u{m)xuIn-M)) • Using a variable p of the coset space, following [151 IS] a-^d [37], we 
express the third equation of p.2p as 



UigMAi) = ('/'A/|f/(5c5«.)l0M)eP''''='^1(^M), g = g^g^., (3.3) 
where we have used the relations 

-^^max 

p=l l<ai<---<ap<M, 
M+l<ii< — <ip<N 

(0a./|C/(5c5^)I0m) = [det(l +pV)]~^ • detw (3.4) 
and the definition 

Piiai ' ' ' Piiap 



•^{Piiai ■ --Pipap) = det 



Pipai ' ■ ■ Pipap 



In (|3.4p a maximum value Mmax is given by Mmax = min(A^ — M, M) and A{- • • ) is an anti- 
symmetrizer. det w; is a determinant of matrix w and is a phase appearing in the decomposition 
of any U{N) matrix as g = g^gw The indices i and a denote unoccupied states (M -|- 1, . . . , N) 
and occupied states (1, . . . , M), respectively. The matrices p and w are defined in Appendix lAl 
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In the Gtm we can introduce an expression cahed the Pliicker coordinate which has played 
important roles for an algebraic construction of a soliton theory in its early stage 



u{9Mm) = Yl •••4lo)> 



l<ai,. ..,aM<N 

9ai,l ■■■ gai,M 



(Pliicker coordinate). (3-5) 



9aM,l ■■■ 9aM,M 

From elementary determinantal calculus, we prove easily the Pliicker coordinate has a relation 

A/+1 

E (-1)'^'4;;:!1m-iA • ^fe;:.'!LiA+i,...,/3M+i = ^ (Plucker relation), 
1=1 

where the indices denote the distinct sets 1 < ai, . . . , om-i < N and 1 < /3i, . . . , /^Af+i < 
The Pliicker relation is equivalent to a bilinear identity equation 

N N 



ciU{gM c^U{gM = ^ U{g)cU(t>) ® U{g)c^\(t>) = 0. 



The bilinear equation has a more general form 

TV N 

Y ciU{gMk) ® CaU{gMi) = ^ U{g)ci\^k) ® U{g)cM = 0, N > k > I > 0, 

where {(pk) and denote A;-particle simple state and /-one, respectively. It is noted that the 
Grj\,f is essentially an SU{N) group manifold since the phase equivalence theorem does hold. 

Now we study the relation between the coset coordinate appeared in (j3.3p and the Pliicker 
coordinates in (j3.5p . Both the well-known coordinates make a crucial role to clarify the algebraic 
relation between the SCFT, i.e. TDHFT, and the soliton theory. 

Using the expressions for unoccupied and occupied states in ()3.3p . we can rewrite ()3.5p as 

ax 

Uig)\(pM) = \M + Yl Y ^t;'^-l,ai+l,...,ap-l,ap+l,...,M,ii,...,v(5C5«;) 

p=l l<ai<---<ap<M, 
M+l<ii< — <ip<N 

^ 4p ■ ■ ■ 4i'^M ■ ■ ■ '^a.p+l^ap-l ' ' ' ' ' ' 

Mn^ax l,...,au...,ap,...,M , X 

I , \ , i,...,M/ X '^i,...,h,...,ip,...,M \yc,aw) 

= I'/'m) +^^i„„,m(5c5«;) 2^ 2^ ^ 

p=l l<ai<-<ap<M, '^l,...,M\9C9w) 
A{[+l<ii<—<ip<N 

X 4iCai---4^CaJ(^M). (3.6) 



The last line of the above is recast again into the form of p.Sp after many time exchanges bet- 
ween • • • Cap and all creation operators so that all the annihilation operators are ordered in 
such a way that they are to the right of all the creation operators including the ones in |(/>m)- 
Then we have the relation 



V 



1,...,M , ^ , ,SrS^-'-''''-'ll,-,au...,ap,...,M 

l,...,ai — l,ai+l,...,ap—l,ap+l,...,M,i 



U-.iMau,) = (-1)^=" ^C^CXZm (9(9: 
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and the following decompositions: 



l,...,ai,...,a 

'h,...,ii,...,ip 

l,...,M 
'h,...,M 



..,M 

,M 



l,...,ai,...,ap,...,M 
h,...,ii,...,ip 



1,...,M/ N N 

^i,...,a/(5'c)^i,...,m(5'»)' 



^i';;;M(5c) = detC(C) = [det(l 



where 



l,...,ai,...,ao,...,M 
■'l,...,il,...,ip 



"m iac) = det 



C(C) 



1,1 



C(C)ai-l,l 
C(C)ai+l,l 

C(C)a,-l,l 

S{C)^p^ 
C(C)ap+l,l 



^^(C)l,A/ 

C(C)ai-l,J\/ 

S{C)ii,M 
C{C)ai+l,M 

C{C)ap-l,M 

S{C)ip,M 
CiC)ap+l,M 



MA 



C{C)m 



MM 



1,...,M 

^i,...,Af vyw 



(g^) = detw. (3.7) 



Here matrix elements in the oi-th, ... and Op-th rows, C(C)ai,i~j\f, • • • and C{C)ap,i,^M are 
replaced with S{C)ii,i'--M, ■ ■ ■ and S{C)ip,i'^M to describe p (1 < p < M) times particle-hole 
excitations from hole state ai to particle state ii,... and those of hole state Up to particle 
state ip, respectively. 

Equating equations (|3.3p and (jS.Sp with equations (j3.6p and (|3.7p , respectively, we obtain the 
anti-symmetrized A{- ■ ■) and the coset variable expressed in terms of Pliicker coordinates as 



l,...,ai ,...,aD,...,M 



•^{Piiai ' ' ' Pipap) 



'1,...A 



(5C) 



1,...,M 
^l,...,Af 



(5C) 



P^a = [S{C)C-\C)]i 



l,...,a,...,M 
h,...,i,...,M 



1,...,M 

h,...,M 



, (3.8) 



in the second Pliicker coordinate of which, only one row matrix elements of its determinantal 
form (j3.7p C{C)a,ir^M are replaced with ^(C)^,!^*/- Expanding the anti-symmetrized A{- • • ) in 
the left-hand side of the first equation of (jS.Sp with respect to, for example, the first column, 
we have a decomposition rule 



l,...,ai ao 



,M 



a,...,ii,...,ip,...,M 



(5C) 



1,...,M 

h....,M 



i9c) 



l)"'"^ Pijai-4(^'iia2 " " ' Pij-i,ajPij+iaj+i ' ' 'Pipap) 

i=i 



^(_l)i+i I'-'^i'-'*^ 



l,...ai,...,M/ N l,...,ai,...,a2,...,aj,...,aj+i,...,ap,...,M , x 
n,\g(^) '^l,„,,ai,...,ii,...,i,_i,...,i, + i,...,v,...,M'.yC^ 



i=i 



l,...,Af 
^l,...,Af 



(5c) 



l,...,Af/ N 

^i,...>/(5c) 



which is rewritten to another form (the second Pliicker relation) 



l,...,Af/ \ l,...,a\,...,ap,...,M f \ 

^i,...,W(5c)^^i,...,ii,...,i,'..,Af (^c) 



, Y^/ l,...,ai,...,A/ / N l,...,ai,...,a2,...,aj,...,aj+i,...,ap,...,A/ . x „ 

+ Z^(-lr^^l,...,i,,...,Af (9c)^l,...,ai,...,n,...,i/_i,...,i,+i,...,V,...,A/(5c) = 0' 



(3.9) 



i=i 



in which hole state oi in the last Pliicker coordinate make no changes (ai — > oi) since in the 
second one particle-hole excitation already occurred from hole state ai to particle state ij [1]. 
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It is weh-known that the Pliicker relation is equivalent to a bilinear identity equation 

N N 

Y,ciU{gMM) CaU{g)\(t>M) = Y,^(9)cl\(t>M)^U{g)ca\(j)M) = 0, 



a=l a=l 

which have made an important role to construct many kinds of solitons on various group mani- 
folds [30]. 

Parallel to the regular representation method by Fukutome \15\ [T6] , we can prove that the 
Lie commutation relation is also satisfied by the differential operators for particle-hole pairs in 
Appendix [B1 

ia A ( * * ^ i9 i * d\ A f ^ ^ i d\ 

d d ^ d d d ^ d d 

From the calculations in Appendix [Bl these differential operators are also proved to satisfy 
relations 

e*"$M,Afb, P*,r) =Pia^M,M{P, P*,t), eai^M,M{P,P*,T) = 0, 

eab^M,M{P,P*,T) = 5ab^Ad,M{p,P* ,t), Cij^ m,M (P^ P* ^ =0' (^-H) 

and a commutator [e*",p*^] = —p*^Pja- A free particle-hole vacuum function (^M,MiP-:P* ^'^) 
given as 

«'M,M(p,P*,r) = [det(l +ptp)]-|e-*^ (3.12) 
Further we can introduce higher order differential operators obeying the relation 

„l,...,ai,...,a^,...,M. o o o N d. i^a^ _ _ _ i^a^, 
^l,...,h,...,i^,...,M P' P ' ' 

which show that by operating the differential operator D on the vacuum function we obtain 
the Pliicker coordinate A. The Pliicker relation (j3.9p becomes a finite set of partial differential 
equations satisfying 

^M,M (p, p* , 'm^^^'A-^^^*' p* > 

p 

+ 2^(-lF^i;...;n,...>^M,M(p,i5 .r)D^^^^^^^^^ .^^^^^^^ ,r) =0, 

l,...,ai,...,ap,...,M / n\* _ / l,...,ai,...,a^,...,M 



(5c5.))* = (^;;:£:::;C:::^(5c)det^ 



'^i,...,ii,...,v,...,M vycy«) 

= ^l,...,ii,...,W..,Af ^M,A/(p,i5 ,r). 

Thus, in both the SCFT and the soliton theory on a group, we can find the common feature that 
the Grassmannian is just identical with the solution space of the bilinear differential equation. 
The solution space of each differential equation becomes an integral surface [32l [3U [1] . 
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3.3 SCF method in 

We will give here a brief sketch of the SCF equation, i.e., the TDHFM. According to Rowe et 
al. [46] . we start with a geometrical aspect of the method in the following way: 

Let us consider the time dependent Schrodinger equation ihdt^ = with a Hamiltonian 

H = hf3ac'pCa + ^{'ja\6f3)c\clci3Ca, (3.13) 

where (7a|5/3) denotes a matrix element of an interaction potential. The starting point for the 
TDHFT lies in an extremal condition of an action integral 

6 r dtC{git)) = 0, C{git)) = {c^M\U{g\t)){ihdt - H)U{gmM- (3.14) 



To get an explicit expression for the TDHFEQ, we calculate an expectation value of one- and 
two-body operators for the S-det (j3.2p . Using the canonical transformation (j3.ip . we have 



M 



Waf3 = {4>M\U{g^)c^f^CaU{g)\(j)Al) = {g'<)f3'p{g^)a'a{4>M\cl,Ca'\4>M) = ^ gaa'gl'f^, (3.15) 

a'=l 

{4'M\U{g'')c\clci3CaU{g)\(l)M) = {g'')y'y{g'')s'5{g'^)f3'l3{g'^)a'a{Hl\Cycl,Ci3'Ca'\4>M) 

= W^^Wps - W^sWp^. (3.16) 

Introducing triangular matrices C{C,) and S'(C) in Gr^ |371ll6j and using an isometric matrix 

nT=[C7T(C),5T(C)], v)u = 1m. 

W in (j3.16p is expressed W = uv) and satisfies VK^ = W (idempotency relation). Then, it 
turns out that the above matrix W is just the density matrix. From p.l6p . we get an energy 
functional, i.e., an expectation value of the Hamiltonian (|3.13p 

H[W] = {(t)M\U{g^)HU{g)\ct>M) = h^aW^p + \ViamW^^Wp5, 

[7a|5/3] = {-ia\5l3) - (7/3|fe). (3.17) 
By projecting the original hamiltonian onto the Grj\//, we obtain also a HF Hamiltonian i?HF[W^] 

HnY[W]=J'^p[W]cicp, jr^p = ^-^^ = h^p + [aP\^5]Wi^. (3.18) 

The Lagrange function C{g{t)) in ()3.14p is computed as 

^{git)) = '-^{g\b9ba + giiiiia - glkdha - g\^g^a) - H[wi (3.19) 

using dtU{g^t))U{g{t)) + U{g^{t)) ■ dtU{g{t)) = 0. The condition (imi) gives the TDHFEQ 

and then we obtain a compact form of the TDHFEQ ihdtg{t) = J-\W {g{t)^]g{t) . The time 
evolution of the S-det ()3.2p is given by 



ihdtU{g{t))\4>M) = HnY[W{g{t))]U{g{t))\(t>M). (3.20) 
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On the other hand, using a u-dependent fermion operator given soon later, from ()3.14p we also 
obtain a compact form of t;-dependent HF equation, instead of time t, as 

ihd^giv) = nW{g{v)}]g{v), {ihd^u{v) = r[W {g{v)]]u{v)) . (3.21) 

Following the observation by D'Ariano and Rasetti's [31] for the relation between "soliton 
equations and coherent states", we may assert that the SCFM presents the theoretical scheme 
for an integrable sub-dynamics on a certain infinite-dimensional fermion Fock space, by iden- 
tifying \(j)M) with the highest weight vector and by regarding the TDHF- manifold Gtm as the 
projection onto a subspace of the r-function. 

We reconstruct a w-dependent SCFM in a Foo and study a relation between soliton equation 
and u-dependent HF equation [JT]. We start from a single-particle Schrodinger equation with 
a u-dependent and a T-periodic potential V{r, v) = V{r, u -|- T) 

hsp{r,v) = -—A + V{r,v), hsp{r,v)'4)a{r,v) = ea'4>air,v). 
2m 

Here we have supposed that an eigen-spectrum is v independent, though the potential is 
dependent on v. It holds an iso-spectrum under a u-evolution of the potential. An eigen- 
function ^(j(r,t;) constitutes an orthonormal complete set and satisfies the same periodicity, 
tjja{r,v -|- T) = Tpaivjv) (Floquet's theorem). This picture is very different from that in [1] 
and [32]. According to Goddard and Olive |47j . we can make Laurent expansion of a fermion- 
field creation-operator ip^{r,v) with a parameter v as 

1 



where z = exp (^271^) given on a unit circle. Thus, the ipNr+a can be regarded as a new fermion 
creation-operator. We obtain also a new fermion annihilation-operator in the same way. The 
anti-commutation relations can be rewritten as 

{ca{v),cl{v')} = dap5iv - v'), {c^{v),cp{v')} = {^(t;), ^(u')} = 0. (3.22) 

Through Laurent expansion of the fermion-field operators, infinite-dimensional fermion opera- 
tors with particle spectra and Laurent spectra can be obtained as 

1 , , 1 



'^(t^-^^') = ^5^exp|i2^^^^^r|, (3.23) 



where Z means the set of the integers. The indices a and r are called the label on particle 
spectra and that on Laurent spectra, respectively. Substitution of (13.231) into ()3.22p leads to the 
anti-commutation relations 



{llj%!r+a,i^Ns+l3} = ^aP^rs, {V'Arr+a- V^W/?} = {V'Afr4«, V^Af^H^g} = 0. (3.24) 

If the canonical transformation (|3.ip has the f-dependence and generates the u-evolution of the 
potential, it is possible to embed a U{N) group induced from (13. ip into a group which can be 
induced from a canonical transformation of the infinite-dimensional fermion operators (j3.24p . 
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According to Kac and Raina |48^ 149] in Appendix [Cl we introduce a Fqo and an associative 
affine Kac-Moody algebra. Here we restrict ourselves to the case of the Lie algebra u{N). The 
corresponding perfect vacuum |Vac) and "simple" state |M) are defined, respectively as 

ipNr+a\yac) = 0, (VaclV'^^+Q, = (r < -1), 
V'TYr+alVac) = 0, {Yac\tpNr+a = (r > 0), 

|M) = V'M---V'i|Vac), (M|M) = 1, (Vac|Vac) = l. (3.25) 

We embed the free vacuum |0) and simple state \(Pm) hito F^o as |0) ^ |Vac), \4>m) > 1-^) 
(M = 1, . . . , A^). Assume that a state with Laurent spectrum corresponding to |0) is a stable 
state with minimal energy. This means a choice of gauge under which |0) corresponds to |Vac). 
The matrix 7 (g u{N)) in U{g) has also the periodicity T. Mapping from a unit circle 
to u{N) [U], we make Laurent expansion of 7 as 7(2) = Yl IrZ^ and impose 7^(-z) = —7(2) 1— > 

7^ = — 7-r and = z* {\z\ = 1). Using the correspondence between basic elements: CaCf^z^ 1— > 
"^{eapir)} = Yl '^N{s-r)+a'4'*Ns+/3 ^-^^ normal-ordered product : 'il^Nr+a'>P*Ns+/3 '■= '^Nr+ai'lfs+p- 
dapSrs {s < 0), let us define the following su(A)(c sl{N)) Lie algebra [Tl [52]: 
X-y = Xry + C ' c, C* = " C (pure imaginary), 

^7 = X] ^Ma,l3 ■ ^N{s-r)+a'<P*Ns\{3 ■, Tr 7^ = 0, 

[X^,c]kM = 0, [X^, XyjKM = ^[7,7'] V) • C, c\M)=1-\M), 

"(7,7) = -a*(7,7') = '^rTT{jrl'-r), (3-26) 

r&Z 

where c denotes a center. As for the t representation (rep) and KM bracket, see Appendix [Cl 
Using equations (|3.24p and (j3.26p . adjoint actions of for ip and i/j* are computed as 

[X^,'4'Nr+a] = 1pN{r-s)+f3ils)l3a, i^-y , i^*Nr+a] = V'^(r.-s)+/3 (7l)/3a ■ (3-27) 

sez sez 

Let us introduce a canonical transformation U{g) = e'^^ satisfying U~^{g) = U{g~^) = U{g^) 
and U{gg') = U{g)U{g'). The g (= e^) has a form analogous to g but with infinite dimension 
and satisfies g^g = gg^ = loo- Further, using (13.270 and the operator identity e^'''Ae~'^^ = 
A+ [X^,A] + ^[X^, [X^, ^]] + • • • , the infinite-dimensional fermion operator is transformed into 

i>Nr+a{g) = U {g)ll)Nr+aU~^ {g) = ^^"^ N{r-s)+l3{9s) (ia, 9Nr+a,Ns+l3 = {gs~r)ap, (3.28) 

where gs means the s-th block matrix of g on each diagonal parallel to the principal diagonal and 
satisfies the ortho-normalization relation. Using the correspondence \4>m) ^ 1^^)) U{g) ^ U{g) 
(= e"^T) and 

N N N 

o=l a=l rSZ o=l reZ 

the bilinear equation on the finite-dimensional Fock space is embedded into the one on F^^ as 

N 

Y,Y.'^Nr+aU{g)\M) (g)rNr-,aU{g)\M) 
a=l reZ 
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N 

= E E UmNr+a\M) Uig)rNr+a\M) = 0, (3.29) 
a=l reZ 

N 

a=l rSZ 

N 

= Y,^U{g)xl^Nr+a\k)'S)Uig)r^,^Jl) =0, k>l, k,l = l,...,N. (3.30) 

where M = 1 ^ N and k > I {k,l = 1 ^ N). Thus we arrive at the following picture: 
An algebra of extracting sub-group orbits made of loop path from Grjvf belongs to the sl{N)- 
reduction of gl{oo) in the sohton theory. Relieving from restrictions of su{N) and (j3.29p and 
taking 7 G sl{N) with M and k > I {e Z), (j3.29p and (|3.30p can be regarded as the bihnear 
equations of the reduced KP (Kadomtsev-Petviashvih) hierarchy and the modified KP in the 
soliton theory [30] . The algebra of extracting various sub-group manifolds made of several loop- 
group pathes [36] from Grjv/ belongs to an s/(A^)-reduction of gl{oo). This picture suggests us 
the possibility to construct a u-dependent HF soliton equation by using the u-dependent SCFM 
governed by conditions (j3.29p and ()3.30p on the space sl{N) bigger than su{N). However we 
must note that in the SCFM bilinear equations (I3.29P and (I3.30p are considered to play roles of 
conditions for assuring of existence of sub-group orbits on Gim different from the soliton theory 
in which boson expressions for them become an infinite set of dynamical equations. Notice also 
that the concept of quasi-particle and vacuum in the SCFM on is connected to the Pliicker 
relation. 

Following m and [32], we embed the original Hamiltonian (13.13P into F^o- By replacing the 
annihilation-creation operators of fermions as ct 1-^ ^ ip^s+is and Cq ^ i^Nr+a^ S^* 

HFoo = h(3a '^Nr+f3i^*Ns+a + ki'y^l^f^) Yl ^Nk+ji^Nl+5^*Ns+f}'>P*Nr+a- (3-31) 

r,seZ k,leZ, r,sGZ 

Previously we had the f-dependent SCF Hamiltonian .7^[VF{<^(f )}] through equations (13.18P 
and (j3.2ip where u-dependence is directly brought from g{v). This is contrast that in the 
usual TDHFEQ (I3.20|) . the time t-dependence of course arises from g{t). To embed this SCF 
Hamiltonian, we introduce a general Hamiltonian on F^o as 

+ ^ E {Nk -\-'y,Nr + a\Nl + S,Ns + /3)ipNk+'Y'^Ni+s'>P*Ns+i3'^*Nr+a, (3-32) 

r,seZ, k,lGZ 

which is equivalent to (|3.3ip if /i7Vs+/3,7Vr+a = and {Nk+j, Nr + a\Nl + S, Ns-'i- P) = {'ya\6l3) 
(equivalence conditions for Hp^) hold. To calculate the formal expectation value of (j3.32p for 
the vector U{g)\M)^ first we do it for one-body and two-body operators. Using (I3.28P we obtain 

(M|V'Ars+/3V'Ar.r+a|-W') = dgr^pa (for r = 0, a = 1, . . . , M and for r < 0, a = 1, . . . , N), 

(for r(s) = 0, a(/3) = 1, . . . , M and for r{s) < 0, a{P) = 1,...,N). 
Then, for one-body and two-body type operators we obtain 

WNr+a,Ns+(3 = {M\U{g^)i^Ms+prNr+aU{g)\M) 
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M N 

= ^{g-r)a'r{9-s)-/l3 + Yi9t~r)a'ri9Ls)ll3^ (3-33) 
7=1 t<0 7=1 

{M\U{g^)ijNk+-y'lpNl+5'llj*Ns+f3'iP*Nr+aUig)\M) 

= WNr+a,Nk+'yWNs+f3,Nl+sWNr+a,Nl+sWNs+f3,Nk+'y 

The W is just the so-called density matrix since it is easily proved to satisfy the idempotency 
relation W'^ = W which has not been given explicitly in [1] and [32]. It provides a strong tool 
to develop our SCF scenario on the -Foo- Taking summation over infinite integers, inevitably 
we have an anomaly in the above expectation value. To avoid this anomaly, the one-body 
operator (I3.33P must be changed to a normal-ordered product as 

{m)ap = {M\U{g^) : r{ep^{-k)}: U{g)\M) = Y.{M\U{g^) : Viv{r+fe)+/J^Wa ^ U{g)\M) 

AI 

= Wi^r+a,N{r+k)+l3 " ^ h,0^l3a = ^ ^(fi'r)a7(5'r„fc)7/3, (3.34) 
reZ r<0 .rGZ7=l 

where we have used the correspondence relation between basic elements. The Wfc is identical 
with a coefficient of the Laurent expansion of the density matrix W ()3.15p 

M 

kez feGZseZ7=i 

We compute formal expectation value of (j3.32p for the state U{g)\M). Introducing a new 
integer K, due to the equivalence conditions for Hp^ , from ()3.3ip we get 

s+a,N{s-K)+l3 

+ h Yl [iK,j),a\iL,6),f3] ^ T?Ar,+ N{r-K)+'y ^ ^Afs+/3,JV(s-L)+5> 
K,Lez r,sez sez 

{hk)f3a = ha, [{k,j),a\{l,6),p] = [ja\6p] y k, I. (3.35) 
Changing W in (j3.35p into its normal-ordered product and using (j3.34|) . we obtain 

kez I lez J 

This result coincides with formal Laurent polynomials of i7[VF] (j3.17p in the sense of the expan- 
sion 

H[W{z)] = J2 I haimU + hba\6p] Y0^i-kUm)/35 \ zK (3.36) 
iez I fcez J 

To get a u-independent Hamiltonian, in (j3.36p it is enough for us to pick up only the term with 
Laurent spectrum / = 0. Then, we may select a density-functional Hamiltonian 

{Hpjm = /l/3a(Wo)a/3 + \ba\5[i] ^ (Wfc)„-, (W_fc);35. (3.37) 

feez 
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This is the extraction of the sub-Hamiltonian -fff^ out of the original Hamiltonian (|3.3ip as 
shown below 

= hf3a^1pNs+f3i^*Ns+a+ lba\^P]Yl Yl '^N{T-K)+'i'^N{s+K)+5'<P*Ns+p'^*Nr+a- (3-38) 

We here adopt ()3.37p as an energy functional for the u{N) HF system on the -Fqo- Through the 
variation 

5{HfJ[W] = Y,{J'-k)o.pKm)po., = Kp5k,o + [a/?l7<5](Wfe)5^, (3.39) 

fcGZ 

we get a SCF Hamiltonian on the F^o similar to formal Laurent expansion of /^hf (|3-18p on the 
GrAf as 

Hf^;BF = ^ ^(-^i<-)a/3 : ll)N{s-K)+a'4^*Ns+l3 ■ ■ (SAO) 

For the f-dependent HF equation on the F^o, the state vector U{g)\M) is required to satisfy 
the variational principle 

6S = 6g dvLig) = 0, L{g) = {M\U{g^){id^ - HpjU{g)\M) , (3.41) 

J Vl 

where we use h = 1 here and hereafter. First by using U{g) = e-^^ we get the following relations: 
5g j dv{M\U{g^)idyUig)\M) = 6g j dv{M\id^\M) + 5g j dv{M\idyX^ 

-^[X^,id^X^] + ---\M), (3.42) 

rez sez 

+ i9„(C-l), (3.43) 

where we have used (j3.28p . From the definition of T{ea/3(r)} and the normal-ordered product, 
we can calculate the ^-differentiation of the second term in the curly bracket of (j3.43p 

: '(/'7v(^_r)+«'0Ars+/3 •= ■ '^{(^ai3{r)} := irdy In z : T{ea^(r)} : . 

sez 

Assume the parameter of Laurent expansion to be z = e~*'^'='". Then (|3.43p is rewritten as 

id^X^ = ^^(A-,;(7r)«/3) : 1pN{s-r)+ai'*Ns+f3 '■ +i9„(C • 1), 

rGZ sGZ 

Dr-v = idv + rujc- (3.44) 

It is seen that in the first term of the r.h.s. of (I3.42P a u-evolution of the reference vacuum 
through z{v) has no influence on variation with respect to g. Concerning the u-evolution 
of Xy (j3.44p a ^-differential acting on ^r{v) and on and ip* through z{v) is transformed to 
a covariant differential D^-v with a connection rcoc which acts only on the 7r(f) from the gauge 
theoretical viewpoint. We denote simply the covariant differential as D. Therefore we can put 
{M\idv\M) = and C = since it has no influence on the energy functional (|3.37p . Then the 
u-differential term in equation (j3.4ip is calculated as 

U{g^)id^U{g) = id^X^ + UidvX^,X^] + X^], X^] + • • • 
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XD^ + Y.U---\^d^Xi^Xii---lXi] + -- - ■ (3-45) 



k>2 



Using (|3.26|) and the symbol D for the covariant differential, each commutator is calculated as 

[id^X^^X^] = X^D-i,'y\ + X] ^ Tr {(i:»7)^7_^} , 
[[idyX^,X^\,X^\ =X[p^^^]^^] + ^rTr{([L»7,7]),.7_,.}, 



[• • • [id^X^, X^l---],X^]= ^] + ^ r Tr {([• • • [D^, 7], • • • ]),,7_,,} . (3.46) 

Substituting p.46p into p.45p and using = D^-v and g = e"^, we get 
U{g^)id^U{g) = X^tDg + C{g^Dg), 

C{g^Dg) = Tr I I ^ • • [i^7,7], • • •],7] | 7-r- i ■ (3.47) 
The expectation value for the reference vacuum is expressed as 

AI N 

{M\Uig^)id^Uig)\M) = ^ ^(gl)„^(D,;„<7s)7a + ^ig^Dg)- 

seZ a=l 7=1 

Using C{g'^Dg) — C{Dg'^ ■ g) = which is proved later, we obtain an explicit expression for the 
L{g) as 

M N 
sGZ a=l 7=1 

Thus L{g) is nothing less than the coefficient of in the Laurent expansion of L{g{z)) (j3.19|) . 
We give another f -dependent HF equation for g: Laurent expansion of 

id^g{v) = J^[W{g{v)]]g{v) and id^U{g(y)m = Hn^iW {g{v))]U {g{v)m . 

Demand the extremal condition of (j3.40p leads to D^g = J-{g)g where J-{g) has an infinite 
A^'-periodic sequence of block form {. . . . . . } like (jC.SP in Appendix [Cl Defining 

{^r)af3{9,^c) =i^cYl s{gsg\-y)a.(5, the D^g = J^{g)g is transformed to 
id^g = J^{g)g, T^ig) ^ Hq) " -^'(5), 



introducing Dy = idv + Hp^.^p, this time which is cast into that on the state vector U{g)\M) as 
D^Uig)\M) = HF^,nFU{g)\M), H'j,^.^^ ^ ^^rU ■■ V'7V(s-r-)+aV'k+/3 
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id^U{g)\M) = H^p^.,,pU{g)\M), H^p^.,,^ ^ ^^?U ■ V'7V(s-r-)+aV'W/3 (3-48) 

r,s6Z 

which suggest symmetry breaking and arising of cohective motion due to recovery of symmet- 
ry. Suppose that g to diagonahze in H^p^.^-p and U(g)\M) to do JT'^ in Hp^.-^-p are 

determined spontaneously when g ~ cpe~'^^'" and dyCp = 0. Using the definition of T'^ we 
have uj(I'{(f') = J'{(p)g^ — (pi where r(^'') has an infinite A^-periodic sequence of block form 
{. . . , -c/° 1, 0,5?, . . . } like ([05]) and 4 = diag{. ..,«,...}. We also obtain grz"" oc g-^^^+'^-'^'v)". 
Thus the quasi-particle energy «(4q,/j = ea^ap) and the boson energy Uc are unified into gauge 
phase. The static f-HFT on Grj\/ has obviously no collective term and leads inevitably to 
u)c^{cp) = 0. (f should compose of only a block-diagonal g^ = , 70 being a block-diago- 
nal su{N) matrix. 

Equation (j3.48p brings a u-evolution of particle degrees of freedom and a common language, 
infinite- dimensional Gioo O'nd affine KM algebra, to discuss the relation between SOFT and 
soliton theory. The SCFT on Fqo is nothing else than the zero-th order of the Laurent expansion 
on Grjv/. Through the construction of the SCFT an explicit algebraic structure of the SCFT 
on Foo is made clear since it is just the gauge theory inherent in the SCFT. The mean-field 
potential degrees of freedom occur from the gauge degrees of freedom of fermions and the 
fermions make pairs among them absorbing a change of gauges. The sub-Hamiltonian (I3.38P 
exhibits such a phenomenon in u{N) algebra, which allows us to interpret absorption of gauge 
as a coherent property of fermion pairs. Thus the SCFTM is regarded as a method to determine 
self-consistently both quasi-particle energy ea{g) and boson energy uJc, to the u-evolution of the 
"fermion gauge". Then we can say that both the energies have been unified into the gauge phase. 

Let e and e* be parameters specifying a continuous deformation of loop path on the Gtm 
and independent on z. Using the notation in (13.260 and calculating in a similar way to ()3.47p . 
^-x^Q^^x^f is obtained as 

^ ^ . ^ C{g-'d,g), 

C{g-^d,g) = a, + 9,7 + E mI- • • [^^7, t], • • • ], t]- (3.49) 

k>2 

To avoid the anomaly, reads 

Y{lr)a(} < Y '^N{s-T)+ai^*Ns+p + ^rO ^ ^a^ > (Tr 7^ = 0). 

rez Uez s<o J 

Then equation (j3.49p is computed to be 

e-^^'d,e^^' = {9r^d,gr)ap ■ V'7V(s-r)+aV'^s+/3 : + Mdo^d^go). (3.50) 

From ^M, we get Cig"^d,g) = Z Tr{g^^d,go) = 0, %^d,go G sl{N,C) and C{g^d,*g) = 

C(ae.gt ■ g) = 0. We obtain also C{g^Dg) = CiDg'' ■ g) = in the same way as the above. 
For subsequent discussion it is convenient to define infinitesimal generators of the collective 
submanifold as follows [501 ttZj '■ 

Xg, = id,U{g) ■ U{gy = Xg, + C{id,g ■ g^) = Xg,, 9^ = id,g ■ g\ 

Xe = id,*U{g)-U{g)^ = Xe + C{id,,g-g^)=Xe. 9 = id,,g-gK (3.51) 
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3.4 SCF method in r-functional space 

Along the soliton theory in the infinite-dimensional fermion Fock space \30\ B9l 154] , we 
transcribe the u-dependent HFT in Foo to the one in r-functional space. We restrict ourselves 
mainly to the cases of sl{N) and su{N) and the group orbit of the fundamental highest weight 
vector \M). Let us consider infinite-dimensional charged fermions ipi and ip* [i € Z) satisfying 
the canonical anti-commutation relation = 5ij and {ip*,Tp1j} = {^/^jj-i/^j} = 0. The 

perfect vacuum |Vac) and the simple state \M) given by (j3.25p are represented in terms of 
another basis Ui (i G Z) and the present fermions ipi and ip* as 

|Vac) ~ z^o /\ ^-1 /\ ^-2 A • • • (a : exterior product), (Vac|Vac) = 1, 
^,|Vac)=0 (?<0), ?/'I|Vac)=0 (i > 0), 
(Vac|V'* = (i<0), (Vac|V'i = (i > 0), 
|M) ~ i/A/ Az^A/^i A--- , {M\M) = 1, 

|M) =V'M---V'i|Vac) (M>0), |M) =VM+i---V'SlVac) (M < 0). 

The basis G Z} is given by the column vector with 1 as the i-th. row and elsewhere. The 
number M is called the charge number. The fermions and ip* [i G Z) generate an algebra 
gl{oo) = {aijipiip*; all but a finite number of a'-s are 0} satisfying j ■(/'*, = ^jk'>Pi'4'i — 
^ii'4'kip*j- We consider further a bigger Lie algebra than gl{oo) so as to include a Heisenberg 
subalgebra (bosons). Following Appendix (Cj it is defined as the vector space Oqo = { X] "^ij 

: tpiip* : +C • c, aij = 0, for \i - j\ > N} and : ipiip* := ipiip* - 6ij ■ c {j < 0). Define a KM 
bracket among such elements Xa = Xa + C • c as 

[Xa, Xb]KM = Xya^b] + a{a, b) ■ c, [Xa, c]km = 0, Xa= ^ aij : ipiip* : . 

For detail see Appendix [O A Heisenberg subalgebra S [49] is defined as 
S = ®kMk + C • c, Afc = ^ : ^i^Pl+k ■■ (keZ). 

From the definition of the normal-ordered product, the boson algebra is obtained as [A^, A/] = 
k6k+ifi ■ c. Afc is called the shift operator and Aq belongs to the center. Suppose level-one 
c = 1. Let be a linear span of semi-infinite monomials with charge number M. For the 

representation of A^ on F^^\ Afc|M) = holds for A; > 0. All the elements A_fc^ • • • A_fcjM) 
(0 < ki < k2 < ■ • ■ < kg) are linear independent with each other. Thus, we have obtained an 
irreducible representation of the algebra 5 in the fermion space F^^\ This is isomorphic to the 
representation of 5 in the corresponding boson space below. 
Let (Tm denote this isomorphism 

ctm: F^^'^ ^B^^'^ = C{xi,X2,...) {deg{xj)=j), |M) ^ 1, 
d 

Ak^—, A^k^kxk {k>0), Ao^M. 

A mapping operator is introduced as um = (M|e^^^\ H{x) = ^ ^j^j {Hamiltonian in r- 

i>i 

FM) [30]. Then we see the correspondence of F^^^) = ©fc>oi^i^''^ with 5^^'^) = ®k&zB[^'^ and 
q(m) ^ z'^'^C{xi,X2,...), where we have defined the direct sum of maps a = (Bm&zO'm and 
have introduced a new variable z to keep a track of the index M and then o"('(/'(^^) = z^ for 
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The contravariant hermitian form on the i?*^*^^ is given as 

where the P* means the complex conjugation of ah the coefficients of the polynomial P and 
X = (XI,X2, . . .). 

We construct a representation in B^^^^ in reduction to sl{N). Let the generating series be 

and introduce Schur polynomials Sk{x) given in Appendix [Dl It should be emphasized that 
in [1] and [32], we already have obtained explicit expressions for the basic elements 



7^ d 

da-. = 



(JM ■ ipiipj '-^ Zij{x,dx), (JM ■■ : ipilp* -.^ Zij{x,dx) {= Zij - 6ij, j < 0) , 
d 1 d 
dxi ' 2 dx2 ' 

Zij{x,dx)= ^ Si+k+fi~M{x)S_j_k+u+M{-x)Sf,{-d.j;)Su{dx), (3.53) 

fj.,u>0, k>0 

which makes a crucial role to construct a f-dependent HFT on U{g)\M) as shown later. For 
any element of the gl{oo) and the Ooo in the B^^'^\ we have got 

au ■ Xa= ^ aijipiil)* ^ ^ aijZij{x,dx), 

i,j&Z 

cjAf : Xa= ^ ttij : ipiip* : +C • 1 ^ aijZij{x, d^) + C ■ 1. (3.54) 

Using exp{i?(x)} = exp { XjA{} = Yl ^jSj{x) and H{x)\M) = due to Aj|M) = 0, 

j>i i> 

x-evolution of the infinite-dimensional fermion operator is given in terms of the Schur polyno- 
mials as 

oo oo 
j=0 j=0 

e^(^V(p)e"^^'^^ = ^(p)e^-^^ e^(^)^*(p)e-^(^) = V*(p)e^-^ (3.55) 

Following Appendixini under the action U{g), the group orbit of the highest weight vector |M) 
is mapped to a space of r-function TM{x,g) = iM\e^'^^^U{g)\M) . Let the Pliicker coordinates 

M, M-l,...,l / N 1 M,M-1,...,1 /X J i I M, M-l,...,l I M,M-1,...,1 j- • i ^ A • * 

t;. . . (g be . a [q] = det o- • '• , Q, , ' ■ : Matrix located on mter- 

section of rows iMjiM-i, ■ ■ ■ and columns M, M — 1, . . . , 1 of 5 and iu = + o etc. The 

Schur-polynomial expression for the r-function is given in a compact form as 

TM{x,g)= 4M^lMlu'tiS9)Sij^t-M,iM (3-56) 

«Af>«A/-i>--->n 

By using (13. 5p and (I3.55p . the derivation of the above is made as follows: 

tm(x, g) = (Vac| • • • V'Me^(-) ^ ^^If^M-Zt. ' ' ' |Vac) 

«A/>«A/-1>--->''1 
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M,M-1,...,1 / X 



u/>«A/-i>--->n 
111 

Inj=M /j,f_i=M 7i=M 

M,M-1,...,1 



^(x) • • • S'ii-/! (x)(Vac|V'r • •V'mV'/ 



•t/'/i |Vac) 
(3.57) 



«A/>«A/-i>--->n 



X det 



«A/ 



-A/ 



(x) 



(x) 



'^«A/-i-(Jl/-l)-l(^) ^iM-i-iM-l)i^) Sij^j_^^(M-l)+l{^) 



ii-h 



-{M-l){x) 'S'ii-l-(A/-l)+l(2;) 'S'ji-l-(A/-l)+2(2;) 



'S'i„-A^+(Af-l)(3;) 
'^iA/-i-(M-l)+(A/-2)(2;) 

'S'ji-1-(A/-1)+{A/-1)2;) 



This equation reads (j3.56p . the generahzation of which to infinite-dimension is given in 

To see that the affine Kac-Moody algebra associated with the Lie algebra gl^ is contained 
as a subalgebra, we give a reduction of gl^ to sin- A subalgebra Xa ( = X] • '^i'^j '■ • l) 

of Ooo is called n-reduced if and only if the following two conditions are satisfied: 

N 

(i) ai+N,j+N = aij and (ii) ^ai^i+Nj = (i G 

i=l 

From (i) and Anj = ^ : V'jV'J'+tvi • 0' ^ [-'^aiATVj] = is proved. This means TM{x,g) 



{g G sl{N)) is independent on xnj, though A^j does not satisfy (ii). As a result, the Hirota's 
equation includes no xjsfj. Adopting the prescription for r-FM, we transcribe the fundamental 
equation (j3.29p for w-dependent HFT on U{g)\Ad) C F^^''^ into the corresponding function 
C B^^^^ in the following forms: 

(1) U{g)\M){U{g) = e^^; G sliN)) ^ TV-reduced KP r-function, 



rM(x,5)7VKP = (M|e^(")C/(5)|M), 



d 



dxN 



■TMix,g)NKP = 0. 



The 5/(00) group symmetry transformation acting on the r-function of the KP hierarchy was 
studied in |5H l3Uj . The KP hierarchy was also studied extensively together with the KdV 
hierarchy by Gelfand and Dickey [52] . The generalized KP hierarchy was obtained out by making 
use of the Gelfand-Dickey approach via the algebra of pseudo-differential operators [53]. The 
r-function TMix, g)NKF bas a f-dependence through g (= e^) in which the anti-hermitian matrix 
7 is given as 7(2;) = ^ 'jrz'' with z = exp (i27rY). 

(2) Quasi-particle and vacuum state 1— > Hirota's bilinear equation (see Appendix [El [301 139]): 



N 



Y,Y.'^Nr+aU{g)\M) (^rNr+aUmM) = 



0=1 



y~] Sj{2y)Sj+N+i{-D) exp j ^ygDs j TM{x,g)NKP ■TM{x,g)NKP = 0. (3.58) 
i>o 



, s>l 



D = {Di,D2, . . . ) denotes the Hirota's bilinear differential operator and D = [Di, \D2, . . . ) 
(3) i;-dependent HF equation on U{g)\M) i— > u-dependent HF equation on rA/(x, ^)7VKp: 



id^U{g{v)]\M) = HFA9{v)}U{g{v)]\M) 
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1-^ idvTM{x,g{v)}NKP = Hf^-uf{x, dx, giv)}TM{x, g{v)} NKP , (3.59) 

in which it is seen that an exphcit and important role of the u-dependence of r-function appears. 
Using (|3.53p and (|3.54p . Hp^-i{-f{x, dx, g) is given as 

HFoo;}lF{x,dx,g) = {J^r{g)}al3ZN(s-r)+a,Ns+(3{x, dx), 

r,s&L 

AI 

{:Fr{g)}ap = Kp5r,o + [a/3|7'5] (W,)^^ , {Wr)aP = Y.^9s)c.^{g\-r)iP ■ (3-60) 

7=1 sez 

3.5 Laurent coefficients of soliton solutions for sl{N) and for su{N) 

We here show typical r-functions called n-soliton solutions. 
On gl{oo) [51]: We get a r-function for gl{oo) as 



rM;n;a,pA^) = (M| 6^^") 6"^^ |M), (3.61) 



which is a famous solution of the KP hierarchy obtained from (j3.58p [54^ HHl H9] . As was shown 
in [21], from the second line of (|3.55p and the Wick's theorem, we get a determinantal formula 
for r-function as 



rM;n;a,p,q{x) = dct < + 



PiJ.-Qu \quj 



If we use T{p, q) t = for a good formal power series of 



TM;n;a,p,q{x) = 6^-^ ^ -1 



V^m/ i<-,3^,<-„ V^m/ \qJ {Pi.- qu){qf,-Pu) 



fj,=l ^ l<fj,<u<n 

where ^,{x,p) = Xjp> and r/^ = C(a;,p^) - i{x,q^) + In (a^tj^r^ 



On Soo: We get a r-function for Uqo as 



On sl{N) [321 HH]: Making a special choice of parameters and (7^ in (j3.6ip as = e*''p^, 
with e = e^'^*/^ and = 1, . . . , — 1, and using [Xa, Aatj] = 0, we get a r-function for sl{N) as 

n 

TAf;n,a,p,e^p(x) = (Mje^^^^e^^i |M). (3.62) 

On 'su{N): Taking the exponent in (j3.62p anti-Hermitian, a r-function for 'su{N) is given as 

TM;n;a,p,.^p(x) = (M|e^(-)e--^ ^ "^^^ ^ |M). 

We have a one-soliton solution on sn(A^) and on the simplest su(2), respectively as 

^{N) : rM;i;a.p,^p(x) = (M|e^(^)e('^^(^)^*(^^^)-'^*<^)'''*(^)}|M) 
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I _|_ g-A/Sg»;(x,p,e''p) _|_ ^Ms ^ri{x,e'' /p* ,1/p*) 



{pp* — l){pp* — 1) ' 
r]ix,p, e'p) = ^XjV (1 - e'^) + Log{a/(l - e")} , 

V (x, e-yp*, 1/p*) = J2 -x,P*-\l - + Log {-a7(l - e"^)} , 
su(2) : TM;i;a,p,-p(x) = (M| e^^^) el'^'^^^^'^* ("^^-'^'H" 5^)^* (^^) } | M) 



2 E xjpJ+Logf -2 E 
1 + e ■'e°'*<^ 



+ 



1)^ 



Finally, we give a reduction of soliton solution in Ooo to the simplest case of sl{2) (KdV). 
A subalgebra Xa (= ^rs '■ '(pri^'t '■ +C • 1) of Ooo and the Chevalley bases for sl{2) [121 EI] 

r,s6Z 

are expressed as 

Xa= (a,.)a/3 : ip2{s-r)+ai^2s+p : +C • 1 (a, /3 = l,2;Tra.r = 0, a.^ G sl{2)), 



r,seZ 

1 




, / 





1 



, h = + 



1 
-1 



(Chevalley bases for sl{2)). 



Using V'(p) = Y i^rP^ and ip{p)* = '^*P~^\ the algebra X" for ?i-soliton in (|3.62p is compu- 

reZ reZ 

ted as 



/ • V'2(s-.)+2pf + /i- • V'2(.-.)+2pf -^)-^VL+2(-P.)^(^^+^)} 



+ 



TV 

E E"^ 

r,s6Z /i=l 



' P^i 



V'2(s-r)+oV'2s+/3- 



0/3 



(3.63) 



Thus, we obtain an n-soliton solution a"(z) = ^ a"^'^' + C • 1 
a matrix as 



alf: = ^^ar{p^,a^), ar{p,a) = a 



_^-2r^ p~(2r+l) 



a^/2) for sl{2) in 



(a^ = const). 



Restricting a solution to the case of Su(2), along the similar way as the above, from p.63p we 
also get 

x:; = x: - xf = x: [a.i^iP.wi-P,) - (-^) r (^) } 

n 

= E )"/3 : ^2(s-r)+a^2s+/3 : + E ("^^^ ~ 
r,sGZ ^l=l 
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which reads 



7r 



We can generahze the above n-sohton solutions to the cases of sl{N) and 'su{N). Using the 
ChevaUey bases for sl{N) and for su{N) \48\ I51|. the Laurent coefficients can be derived for 
each case as 



soliton 



J2 (Pm' e^^^M' ^t^)+^Y 



sl{N) : <| 
a,{= a, {p, e-'p, h)^^ (a, /? = 1, . . . , iV)} = 



~ hp 



1 



P 



P 



-(7V-1) 



Afr+a-/3g-s/3\ 



-2s 



~Ns 



7r{= 7r (P, e'^P, (a, /? = 1, . . . , iV)} = flr 



*T 



- e 1 - e* 



-Nj- 



*Nr 



6* J V 7 
P 



1 
P 

N-1 



P 



P 



-2s 



p^ 



p- 



-(^-1) 



p 



*(iV-l) 



P^ 



P 
1 



Suppose the external parameter f to be a time t and solve the TDHFEQ by restricting 
a solution space to the above spaces. We get a soliton solution, i.e., a solitary wave propagating 
on a surface rather than a colliding soliton [55]. This is in contrast with a two-dimensional 
soliton [56], i.e., dromion [S^, of the Davey-Stewartson equation (DSE) [58] which provides 
a two-dimensional generalization of the celebrated nonlinear Schrodinger equation (NLSE) |57j . 
The dromion for gl(2oo) was derived from the standpoint of a Clifford algebra with generators 
of infinite free fermions, in terms of a reduction of the two-component KP hierarchy [59 ^ 160 ^ 161]. 
As Kac and van der Leur pointed out [61], the dromion solution of the DSE was first studied 
from the point of view of the spinor formalism by Heredero et al. [62] . 



3.6 Summary and discussions 

We have transcribed a bilinear equation for i;-HFT into the corresponding r-function using 
regular representations for groups ^ and Schur polynomials. The concept of quasi-particle and 
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vacuum in SCFT is connected with bilinear differential equations. So far SCFM has focused 
mainly on construction of various types of boson expansions for quantum fluctuations of mean- 
field(MF) rather than taking the bilinear differential equations (Pliicker relations) into account. 
These methods turn out to be essentially equivalent with each other. Various subgroup-manifolds 
consisting of several loop-group paths [36] exist innumerably in Gr^ relating to collective motions. 
To go beyond the perturbative method in terms of the collective variables, we have aimed to 
construct f-HFT on affine Kac-Moody algebras along soliton theory, using infinite-dimensional 
fermions. These fermions have been introduced through Laurent expansion of finite-dimensional 
fermions with respect to the degrees of freedom of fermions related to MF. Consequently u-SCFT 
on Foo leads to dynamics on infinite-dimensional Grassmannian Groo- Grjv/ is identified with 
Groo affiliated with a manifold obtained by reduction of g'^oo to sl{N) and sii( A^) (reduction of 
KP hierarchy to DS, NLS and KdV hierarchies). We have given explicit expressions for Laurent 
coefficients of soliton solutions for sl{N) and 'su{N) using Chevalley bases for sl{N) and su{N). 
In this sence the algebraic treatment of extracting subgroup-orbits with z{\z\ = 1) from Gtm 
exactly forms the differential equation(Hirota's bilinear equation). The u-SCFT on Fqo results 
in gauge theory of fermions and collective motion due to quantal fluctuations of u-dependent 
SCMF potential is attributed to motion of the gauge of fermions in which common gauge factor 
causes interference among fermions. The concept of particle and collective motions is regarded 
as the compatible condition for particle and collective modes. The collective variables may have 
close relation with a spectral parameter in soliton theory. These show that u-SCFT on F^o 
presents us new algebraic method on for microscopic understanding of fermion many-body 
systems. 

We have studied the relation between f-SCFT and soliton theory on group manifold and 
shown that both the theories describe dynamics on each Grassmannian Gr which is the group 
orbit of highest weight vector. The former stands on the finite-dimensional fermion operators 
but the latter does on the infinite-dimensional ones. Each Gr is just identical with the solution 
space for respective finite and infinite set of bilinear differential equations on the boson space 
mapped from those on the fermion space. We have investigated the dynamics on f-dependent HF 
manifold using regular representations for groups [16]. A picture of quasi particle and vacuum 
in u-SCFT is connected with the bilinear differential equation. u-HFT on finite-dimensional 
Fock space is embedded into u-HFT on infinite one. The wave function in an SC T-periodic MF 
potential becomes dependent on the Laurent parameter z on a unit circle S^. This owes to the 
introduction of affine Kac-Moody algebra by infinite-dimensional fermion operators, Laurent 
expansion of finite-dimensional fermions with respect to z. The Pliicker relation on coset variab- 
les becomes analogous to Hirota's bilinear form. The u-SCFM has been mainly devoted to the 
construction of boson-coordinate systems rather than that of soliton solution by r-FM. It turns 
out that both the methods are equivalent with each other due to the Pliicker relation defining Gr. 
From loop group viewpoint and with clearer physical picture we have proposed description of 
particle and collective motions in f-SCFT on Fqo rclcitioii to iso-spGctro^l eQuation in soliton 
theory. Then the i;-SCFT on F(y^ may be regarded as soliton theory in the sense that it bases 
on Groo and may describe dynamics on infinite set of real fermion-harmonic oscillators though 
the soliton theory describes dynamics on complex ones. The soliton equation is nothing but the 
bilinear equation and the boson coordinate Xk with highest degree plays a role of an evolutional 
variable on r-functional space (FS) on which in f-HFT, the bilinear equation provides algebraic 
means to extract subgroup orbits parametrized with z from Gr^f . The infinite set of Xk becomes 
coordinates on r-FS and their u-evolution yield trajectories of the SCF Hamiltonian Hp . 

Though we have started with a periodic potential to introduce infinite-dimensional fermions, 
it is easy to see that f-dependence with periodicity T is by no means a necessary condition. The 
fact that Schrodinger function is dependent on an unit circle S^, however, makes a crucial role 
for construction of infinite-dimensional fermions. As pointed out in [32], it turns out that the 
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fully parametrized u-dependent SCF Hamiltonian is made up of only the i;-dependent Hamilto- 
nian Hp^-np. Then, we have a very important question why infinite-dimensional Lie algebras 
work well in fermion systems. As concerns this problem, Pan and Draayer (PD) [63] have de- 
veloped an infinite-dimensional algebraic approach using affine Lie algebras su(2) and su{l, 1). 
They have introduced fermion pair operators with two parameters for the general pairing Hamil- 
tonian and boson operators through Jordan-Schwinger fermion-boson mapping for an exactly 
solvable su{2) Lipkin-Meshkov-Glick (LMG) model [32]. They have obtained analytical ex- 
pressions for exact eigenvalues and eigenfunctions of this Hamiltonian based on the Bethe an- 
zatz (BA), from which BA equation [63] or Richardson equation |65j is derived. 

It is interesting to study a relationship between various subgroup-manifolds of Groo and col- 
lective sub-manifolds of u-SCF Hamiltonian by using a simple and exactly solvable LMG model. 
Notwithstanding, it is possible to provide a theoretical frame of formal RPA [34^ [35] as a tool 
of truncating a collective motion with only one normal mode, i.e., a collective submanifold out 
of Groo. As mentioned in [341 135] . the collective submanifold may be interpreted as a rota- 
tor on curved surface in Groo- It is stressed that the u-HFT on Foo describes a dynamics on 
real fermion-harmonic oscillators while soliton theory does the same but on complex oscilla- 
tors. This remark gives us an attractive task to extend the u-HFT on real space 'su{N) to the 
theory on complex space sl{N,C) removing the restriction \z\ = 1. We have discussed a close 
connection between u-SCFM and r-FM on an abstract fermion Fock space and denoted them 
independently on S^. It means that algebro-geometric structures of zn/inite-dimensional fermion 
many-body systems is also realisable in yinite-dimensional ones. The u-dependent HF equation 
on TM{x,g), however, should lead to multi-circles, relating closely to a problem of construction 
of multi-dimensional soliton theory [661 [67l [61] . It is also a very exciting problem to investi- 
gate such new motions on the multi-circles (d: Number of circles) in finite fermion many-body 
systems. As suggested by the referee, the motions, on the other hand, may be related to the 
coupled (d + 1)D systems [61] or the linear flows on the Birkhoff strata of the universal Sato 
Grassmannian [68] . 

4 RPA equation embedded into infinite-dimensional Fock space 

4.1 Introduction 

The purpose of this section is to give a geometrical aspect of RPA equation (RPAEQ) [50\ [69] 
and an explicit expression for the RPAEQ with a normal mode on i^oo- We also argue about 
the relation between a loop collective path and a formal RPAEQ (FRPAEQ). Consequently, it 
can be proved that the usual perturbative method with respect to periodic collective variables rj 
and if in TDHFT [25], is involved in the present method which aims for constructing TDHFT 
on the affine KM algebra. It turns out that the collective submanifold is exactly a rotator 
on a curved surface in the Groo- If we could arrive successfully at our final goal of clarifying 
relation between the SCFT and the soliton theory on a group, the present work may give us 
important clues for description of large-amplitude collective motions in nuclei and molecules 
and for construction of multi-dimensional soliton equations |61[ [66] since the collective motions 
usually occur in multi-dimensional loop space. 

4.2 Construction of formal RPA equation on F^o 

We construct the FRPAEQ on F^o- We put the following canonicity conditions which guarantee 
the variables (e,e*) to be an orthogonally canonical coordinate system [6l [25l [32] : 



m\g) ^ {M\U{g^)d,U{g)\M) = \e\ {g\d,*\g) ^ {M\U{g^)d,,U{9)\M) = -h- (4-1) 
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Previously we define the infinitesimal generators of the collective submanifold and Xq^ 
(|3.51|) in which the term C{g~'^d^g) is proved to vanish. From these infinitesimal generators and 
de* {g\de\g) — de{g\d^*\g) , we obtain the weak orthogonality condition 

M N 

1 = {g\[Xe,Xsm = ^ J]5]([0,0t]^.)a7(W-r)7" + E^Tr(Ml.), (4.2) 

a=l7=lrGZ rSZ 

where we have used (|3.26p and (|3.34|) . 

As shown in [32], using Lax's ideas [70] we recast (|3.47p and Dtg = J-{g)g, and (j3.5ip into 

Dtg = T{g)g, dtg'' = 0, J^{g) = Hf), 
ideg = 9\g)g, 0\g) = e\g^) + g\d,i)f^ • t, 

id,*g = 0{g) = e{f) + g\dy,)f^ ■ t. (4.3) 

M 

Upon introduction of E = ^ eQ(e,e*), the canonicity condition (|4.ip transforms into 
{9\de\9) = - id,E ■ t = ^e* - id,E ■ t, 

{9\deA9) = {f\de*\f)-id,,E-t = -\e-id,*E-t. (4.4) 
From (j4.4p . the weak orthogonality condition (j4.2p is expressed as 

1 = d,*{m9)-dMdA9) = = (/|[X,(g«),X,t(go)]|5°). 

To satisfy integrability conditions for e, e* and i, curvatures obtained from (j4.3p should vanish; 

Ct,, = Dte\g) - id^Hg) + [0\g),:F{g)] = o, 
Ct,,. ^ Dte{g) - id^.Hg) + m),Hg)] = o, 

C,,,. ^ id,e{g) - id,,e\g) + [e{g),e\g)] = 0, (4.5) 

and dtg^ = 0. Here DfO and DtO"^ are defined as 

{Dte)r = Dr.tOr = {idt + rLO,)er, {DtO^^ = Dr.tOl, = [iSt + ruo,)el,. 

The expressions for the curvatures on the quasi-particle frame (QPF) have the same form 
as those of RPAEQs in the finite Fock space [5^. As mentioned before, the TDHFEQ on 
the leads to the RPAEQ if we take into account only a small fluctuation around a stationary 
ground-state solution. The form of RPAEQ on the QPF has a following simple geometrical 
interpretation: Relative vector fields made of the SCF Hamiltonian around each point on loop 
pathes also take the form of RPAEQ around the same point which is in turn a fixed point in 
the QPF. Thus, the curvature equation in the QPF is regarded as the FRPAEQ on the Grgo- 
Using (I3.28p . the canonical transformation for g is given by 

TV 

seZ/3=i 

together with its hermitian conjugate. Owing to [50], ()4.3p is rewritten on the above QPF as 

-Dtg^ = J-(5^)|qpf5^ -^(5^)iqpf = g'^nm, 
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(^"^)lqpf5^ 



-id.^g^ = eig%pig\ Oig^ 



Iqpf = 9 



(4.6) 



The subscript "qpf" means the quasi-particle frame (QPF). For (j4.5p we obtain also another 
expression on this QPF as 



{DtO^ - ideT-[e\T])\^pf = 0, (DtO - id,*J'-[9,J'])\^pf = 0, 



(4.7) 



Further, using (j4.6p and the relation ideJ-'\qpi = 'ide{g'^ T{g)g) = — [0^, .7^]|qpf + gUde^g, one can 
rewrite equations in the first line of (j4.7p as 



A^^lqpf - g^id^nm = 0, A^lqpf - gUd,,J'ig)g = 0. 
From (j4.6p and (j4.3p . the infinitesimal operators are expressed as 
^^(5^)|qpf = -id J • g = e'''{dA • t + et(^ot)|^p^|e-.a^ 



(4i 



(4.9) 



together with the same relation for 0(^"^)|qpf. We have also 9^{g^^)\qp{ = —id^cp^ ■ ip and 
6{cp'^)\cg,f = —ide*g^^ • g^. Then, from (j4.8p we can derive the FRPAEQ on the Groo in the form 



{^^(^°^)lqpf} 



^,^^(5°^)lqpf 



^9 



^^d,nf)9' = 0, 



(4.10) 



To obtain an explicit expression for the last term of the l.h.s. of (I4.10p . we introduce an auxi- 



liary density matrix R = ^'^diag [• • • / 



M(Si{N-M) 



g^\ where / 



M(SiiN-M) 



-hi 



LN-M 



The R is related to density matrix as ii = / — 2W {I: infinite-dimensional unit matrix). 
Then, we obtain 



- id.g^'' ■ (flM®{N-M) - hm{N-M)i-ideg°'' ■ 5°)}5°'^ 
i9'[0\9'%pi,hmN-M)]9'^, 



id,W 



and we have used (j4.9p . Further we introduce the following quantities: 



(4.11) 



qpf 



bI 



qpf 



d 1 



r''^lqpf> Im®{N-M)\ 



-00 



Using these, we rewrite (j4.1ip as 

id^Wr = J2 9lBU_XpigV = E 



9l 



k,l£Z 



k,iez 







9i 



Pt 



(4.12) 



Let a (a) and i (i) be 1, . . . ,m hole-states and m + 1, . . . , particle-states of the QPF, 
respectively. Substituting the second equation of (|4.12p into (|3.39p . for r / we get 



k,iez 
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Thus, we can reach the desired form of the equation, part of the FRPAEQ on the Groo (|4.10p . 



E 





E 


gV 








kl 




kl 


1 F\ 




(xh 




ici 


kl 


\D\ 


kl 




ia 




ia 


kl 


\D\ 


kl 




ai 




ia 


kl 




kl 


\F\ 








ia 



^{k-l)-{k-l)-r 
^{k-l)-{k-l)-r 



^(k-l)-{k-l)-r 
^{k-l)-(k-l)-r 



kl 




kl 




\F\ 




ah 




ai 


kl 




kl 




\D\ 




ia 




ai 


kl 




kl 




\D\ 




ai 




ai 


kl 




kl 


\F\ 




ij 




ai 



^{k~l)~{k~l)-r 
\k-T)-{k-l)-r 



6°- - 
{k-l)-{k-l)-r 



6°- - 
{k~l)~(k~l)~r 



Substituting the above result into (j4.10p , we can derive the FRPAEQ on Fqo ■ 

Finally we show the following equations to determine the collective submanifold and motion: 
The canonicity condition ()4.ip : 



M 



a=l s& \ 



The FRPAEQ (ITOD : 



0, 



5 = 5°(e,e*)e-«(^'^*)*. 



Through constructions of the TDHFT and the FRPAEQ on Fqo, the following become apparent: 
The ordinary perturbative method for collective variables t] and rf [25] is involved in the way of 
construction of the TDHFT on the affine KM algebra if we restrict ourselves to 'su{N). When 
the rj and r/* are represented as rj = we can always express 



E 



on the Lie algebra if we put z = e^'^. This means that the infinite-dimensional Lie algebra 
in the SCFT is introduced in a natural way and is useful to study various motions of fermion 
many-body systems. 



4.3 Summary and discussions 

FRPAEQ has been provided as a tool for truncating a collective submanifold with only one 
normal mode out of an Groo- We have given a simple geometrical interpretation for FRPAEQ. 
The collective submanifold is interpreted as a rotator on a curved surface in the Groo- In F^o, 
to study motions of finite fermion systems, it is manifestly natural and useful to introduce an 
infinite-dimensional Lie algebra arising from anti-commutation relations among fermions. In 
order to discuss the relation between TDHFT and soliton theory, we have given expressions for 
TDHFT on r-FS along soliton theory. From the loop group viewpoint and with a clearer physical 
picture, we have proposed a way of describing particle and collective motions in SCFT on F^o 
in relation to an iso-spectral equation in soliton theory. Then, SCFT on Fqo may be regarded 
as soliton theory in the sense that it is based on the Groo and may describe dynamics on an 
infinite set of real fermion-harmonic oscillators. On the other hand, soliton theory describes 
dynamics on complex fermion-harmonic oscillators. It is one of the most challenging problem 
to extend real space su{N) to complex space sl{N) in TDHFT on Fqo tog6th.Gr with roniovcLl 
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of the restriction \z\ = 1. Concerning the construction of sohton theory on multi-dimensional 
space |6H I66j . we have an interesting future problem that is to extend the Pliicker relation 
(Hirota's form) with only one circle to the case of multi-circles such that SCFM on F^o can 
describe dynamics of fermion systems in terms of multi-RPA bosons. 

5 Infinite-dimensional KM algebraic approach to LMG model 

5.1 Introduction 

To go beyond the maximaly-decoupled method, we have aimed to construct an SCF theory, 
i.e., u-HFT. In constructing the v-HF theory we must observe, however, the following two 
different points between the maximaly-decoupled method and the t;-HF SCFM (i) The former 
is built on the finite-dimensional Lie algebra but the latter on the infinite-dimensional one. 
(ii) The former has an SCF Hamiltonian consisting of a fermion one-body operator, which 
is derived from a functional derivative of an expectaion value of a fermion Hamiltonian by 
a ground-state wave function. The latter has a fermion Hamiltonian with a one-body type 
operator brought artificially as an operator which maps states on a fermion Fock space into 
corresponding ones on a r-FS. Toward such an ultimate goal, the f-HFT has been reconstructed 
on an affine KM algebra along the soliton theory, using infinite-dimensional fermion. An infinite- 
dimensional fermion operator is introduced through a Laurent expansion of finite-dimensional 
fermion operators with respect to degrees of freedom of the fermions related to a u-dependent 
potential with a T-periodicity. A bilinear equation for the u-HFT has been transcribed onto 
the corresponding r-function using the regular representation for the group and the Schur- 
polynomials. The f-HF SCFM on an infinite-dimensional Fock space -Fqo leads to a dynamics 
on an infinite-dimensional Grassmannian Groo and may describe more precisely such a dynamics 
on the group manifold. A finite-dimensional Grassmannian is identified with a Groo which is 
affiliated with the group manifold obtained by reducting (7/(00) to sl{N) and su{N). We have 
given explicit expressions for Laurent coefficients of soliton solutions for sl{N) and 'su{N) on 
the Groo using Chevalley bases for sl{N) and su{N). As an illustration we make the u-HFT 
approach to an infinite-dimensional matrix model extended from the finite-dimensional su{2) 
LMG model and represent an infinite-dimensional matrix LMG model in terms of the Schur 
polynomials. 

5.2 Application to Lipkin— Meshkov— Glick model 

To show the usefulness of the infinite KM algebra and to avoid an unnecessary complication, we 
apply it to a simple model, the LMG model consisting of N (= M) particles. Let us introduce the 
LMG Hamiltonian which has two A^-fold degenerate levels with energies and — ^e, respectively 



where the indices i and a stand for particle-state and hole-state, respectively and satisfy the 
SU{2) quasi-spin algebra 



H = eKo-lV{Kl+K^^). 



The operators Kq, and A'_ are defined by 




(5.1) 



[K+,K_]=2Ko. 
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Then the S-det, jS*^), in which the N particles fill the lower level, satisfies 



Ci\S 



0, 







{i,a = l,2,...,N), 



K\S 



0. 



We here use a notation [c^] denoting a 2A^-dimensional row vector [cji,c|] {i = 1,2, ... ,N] 
a = 1,2, ... , N). We introduce the following SU{2N) Thouless transformation: 



Uig)[c^]U-Hg) 



9 



cos I • e~^5('^+v') -1^ - sin f • e'^^C^-^) • Ijy 



sm ^ ■ e°2'-'^ ■ In 



l2N, det5r=l, 
Uig)Uig') = U{gg'), U{g~') = U-\g) = U\g), 



9^9 = 99^^ 



(5.2) 



U{1) 



1. 



The above SU {2N) matrix is essentially the direct sum of the SU{2) matrix. Any A^-particle 
S-det is constructed by the Thouless transformation of a reference S-det, \S^) (the Thouless 
theorem) as 



g) = U{g)\S^) = {S''\U{g)\S'')e^v\p^Mim+]\S 



(5^|C/(<7)|5^) = (cos f f , p= {p,a) = tan (f ) e'^ ■ In, 

which is the CS rep of fermion state vector on the SU{2N) group [5]. 
The HF density matrix is given as 



(5.3) 



W 



COS^ I • 1 AT 



■ sm c/e 



N 



IsmOe'^f'-lN sin^f-l^v 



= ^l2N + cose\h2N + ^smOe ''''e2N + ^sinee'^^hN, 

where I2N is a 2A^-dimensional unit matrix and h2N, (^2N and f2N are defined in the next 
subsection. The usual HF energy {g\H\g) (= if[M^]) is obtained as 



H[W] 



eN 
2 



sin^l -cos^l -x{^sin0e-^^}^ -Xl^sin^e*'^} 



X 



(jv-i)y 



The Fock operator (= 5//[VF]/5VF^), in the HF approximation is represented as 



F[W\ 



sm 



6le-# . 1 



N 



-6X2 ^™ ■ 

■^-Ij 



£7 



■^h2N 



f sin Oe'^f'siN - f sin Oe-'^fiN- 



(5.4) 



5.3 Infinite-dimensional Lipkin— Meshkov— Glick model 

Using (I3.25p . we introduce the following infinite-dimensional "particle"- and "hole" -annihilation 
operators V'ATr+j particle-state and i^Nr+a for hole-state, respectively, and a vacuum sta- 
te \M): 

rNr+i\M)=0 {i = l,...,N, r>0), ^PNr+a\M) = ia = l,...,N, r<0), 
|M) = VM---^i|Vac) {M = N). 

Then the Kq and K± (15. Ih are expressed in terms of the above infinite-dimensional operators 
as follows: 

^0 = ^ (^/l2Jv) : '(pN{s-r)+a'4'*Ns+l3 ' 
r,s€Z 
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/ N N \ 

r,seZ \j=l a=l / 

AT 

= X^ (e2Ar)a/3 : '4'N{s-r)+atp*Ns+/3 X^ X^ ' (s-r)+j V^A^s+a ^ 
r,seZ r,s£Zi=a=l 

N 

-f^- = X^ {f2N)af3 ■ 'ipN{s-r)+a'^*Ns+l3 •= X^ ' '^Af(^*-r)+a'0Ars+j 



r,s£Z a=i=l 



and 



TV 



^_|M) = 5^ ^ : i^Nis-r)+arNs+i ■ \M) =0 (s > 0). 



r,s£Z a=i=l 



Let us introduce the following 2A'^-dimensional dual elements of the direct sum of the algebra 
sl{2,C) multiplied by z^': 



Ko{r) = IhiNz" = \ 



K+{r) = e2NZ 
(r) = f2NZ 



In 
-1^ 

1^ " 



In 



= e(r), 
= fir); 

Ko{r) = {ko{r)} ^^eafs = (^fh2N z""^ 6^/3 = (^W) e„/3(r), 



f e2N^ 



2N^ J af3 



/ e2N 

\f2N/ af3 



e«/3(r). 



Using the formulas in Appendix[Cl the T reps of the operators Kq and K± are given, respectively, 
in the following forms: 



2h2N) r{eap{r)}, 

/ ap 



2N I EN(^g_j.)+a,Ns+l3 



a/3 



/J (^^27v) ^'^N{s-r)+a'^*Ns+l3i 

^ ^ ^ ' OLp 



re2 



\]2N 



a/3 



E 



/ e-2N 



E 



2NJ a[3 



N{s-r)+a,Ns+i3 



E 

r,seZ 



(^2n\ I 

? VN(s-r)+aVNs+B^ 
J2NJ afi 



from which we have the KM brackets among the operators K^ir^ and K±{r'). For detailed 
calculations see Appendix [Fj 

[i^o(r),i^o(s)]KM = h^rbr^sfi ■ c, [K±{r),K±{s)]^y^ = 0, 

[Ko{r),K±{s)]^^ = ±K±{r + s), [K+{r), K^{s)]^^ = 2K^{r + s) + Nr6r+s,o • c. 

Following Lepowsky and Wilson [7l], we introduce the elements of s/(2,C) in the linear- 
combination forms of e2N and f2N- 



= {e2N)al3eaf3{r) ± {f2N)al3eal3{r + 1) 



Yk 



a/3 
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In ■ 









In ' 


_ zIn 







. -zIn 








= 

the 2-cocycle a's on a pair of the above elements read 



z'' (/c = 2r + l), 



OL 



a 



N 



N 



^ Xk,-^Yi 
N ./N 



N 



Sk+i,o, a ( 



—kSk+ifi, 
-^k+i,o 



{k = 2r + 1, I = 2s + 1). 
For details see Appendix [Fj Then we have the following KM brackets and the map ax- 



Xk, —?^Xi 



N 



N 



^ Xk,4^Yi 



N VN 
1 



ax : 



ax ■■ 



Xk 



KM 



KM 

d 
dxk 



kxi 



^ Yk,^Yi 



N 



N 



-kSk+ifi ■ c, 



KM 



1 



NVN 



Yk+i + 4+i,o • c, 



1 



N 



^^Yk ^ kyk, 



N 



d 

dyk 



Then ^^KQ{r), -^Xk {k = 2r + 1) and -^Yk {k = 2r + 1) are clearly an infinite-dimensional 

Heisenberg subalgebra of the KM algebra sl{2,C). We also introduce the element in the form 
of h2N as 



Y2r = {i-h2N)af3eal3ir)) = {Y2r)al3eaf3, 5^2r 



-In 
In 



5.4 Representation of infinite-dimensional LMG model 
in terms of Schur polynomials 

The expressions for the operators Kq and K± (jS.ip in terms of the operators -^X2r+i, -^^2r+i 
and -^Y2r and the expressions for the map ax for the operators -^X2r+i and -^X_(^2r+i) 
are given as follows: 
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E 
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X- 



N 

K+{r) 
K_(r + 1) 
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N 
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2r+l 



aK ■■ 



X_ 



(2r+l) 



(2r + l)x2r+i (r>0). (5.5) 



Consequently we can obtain important sum-rules for the operators l2r+i and Y2r as 

Y2r _ 

J ~ ^[So^ , 1 (2x) 1 

s>0 

-2 Yl ^2m. + l 2 Yl 277+1 SiOn+l 
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(5.6) 
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r+l 



4E ^2^+1 -4E 2^ 



+1 ax. 



4E a:2m+i 4E aSTT 



+ 1 dx 



2n + l 



(5.7) 



where we have used the relation 



(X2n = 0, n > 1), 



have an expression for a quadratic operator + as 



which is derived from the definition of the Schur polynomial Sk{x) in Appendix [Dl Further we 

2 r \2 

1 \ I 

--X2r+l + - 



Finally from (|5.5p . (j5.6p . (j5.7p and (jS.Sp we get an expression for the LMG Hamiltonian as 
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-^|exp(4X;.2...)exp^-4g^^j 
-exp(-4|;.2...)exp(4|^^)}. 



(5.9) 



5.5 Infinite-dimensional representation of SU{2N)oo transformation 



We prepare operators KQ{ip), Kq{'iIj) and K±{9) to generate an infinite-dimensional representa- 
tion of an SU{2N)oo transformation-matrix. First we give Kq{lp) and Kq{iP) as 



-<^r \ 1 
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adjoint actions of which for -i/'Afr+a and ipp^^+a computed as 



f '4'Nr-^\ 



2^2iV 



/3q 



Then the infinite-dimensional fermion 



operator 7pNr+a{g(^'p-^) 



is transformed into 



E 



fi^N(r-s)+f}{9(Vs-^^) 



9(Vs 

\i>3 



e y-ipsJ 2 



(5.10) 



In a similar way as the above we also give K±{6) as 

K+{0) = k^_{e) = ^ {-0re2N)o.p ■■ i'N{s-r)+ai^*Ns+l3 ■ 




1^ 



: 'lpN{s-r)+a'^*Ns+j '■ '■ N (s-r)+a'^*N s+b 



N 



= E E • '^N{s-r)+ii'*Ns+a 
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and from (|I.2p adjoint actions of which for ipj^r+a are given as 



K-iO) 



,'4'Nr- 



s-r)+/3 



-0-rf2N/ f^o 



{Or = e_r). 



Then the transformed infinite-dimensional fermion operator ipNr+a{ge) is transformed into 



E V'Af(r-s)+/3(5£'s)/3a, gOs = 



hie2N-f2N)i-es) 



(5.11) 



The whole transformation of the infinite-dimensional fermion operator ^Nr+a{g) is given as 
i'Nr+a{g) = U{g^)U{g0)U{g^)i^Nr+aU~^{g^)U~^{ge)U~^{g^) 

= U{g^)^U{gg)ll)Mir-s)+f5U^^{ge)U^^{g^) {g^s ) (Sa 



E u{g^)TpN{r-s-t)+-yU ^{g^){get)-ip{gvs)po 

E '^N{r-s-t-u)+S {g^pu )S'r (9et)'rl3 {g^s )i3ol 
E '^N(r~s~t~u)+I3{gu,t,s)l3a- 



s,t,u 



(5.12) 



S,t,U(: 



The block matrix gu,t,s of an SU{2N)oo transformation-matrix is given by (jL5|). For simplicity, 
in the sum over u and t of (I5.12p . we pick up only the term u = s and t = —s, then we have 
a simple transformation of the infinite-dimensional fermion operator ipNr+aig) as 

IpNr+aig) = y^^'ipN{r-s)+f3{gs)(3a, 
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5.6 Representation of infinite-dimensional HF Hamiltonian 
in terms of Schur polynomials 



Using the t rep of (j5.4p and (j5.5p . an infinite-dimensional HF Hamiltonian for the LMG model 
is represented in terms of the Schur polynomials as 

= -Yl [^U^^N)ap + £X^smee'^{e2N)ap + ex^smee~'^{f2N)ap T{e„^(r)} 
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the second line of which has a form very similar to the Hamiltonian given by Mansfield |72j . 
The HF density matrix for the LMG model is also represented in terms of KQ{r) and K±{r) as 

W - i/27V = [cos6li(/i2jv)a/3 + ^sin6'e"*'^(e27v)„/3 + ^ sin6'e*'^(/2Ar)a/3 r{eaf3{r)} 
= ^cosOKoir) + lsm9e-'^K+{r) + i sin0e*^ir_(r) 



The infinite-dimensional HF operator for the LMG model Hp^.iip{x,dx,g), corresponding to 
(j3.60p . is expressed in terms of Z2N{s-r)+a,2Ns+i3{Xy dx) {a, P = 1,2, . . . , N, N + 1, . . . , 2N) as 

HF^;UF{x,dx,g) = V |-F[VF](r)| Z2N{s-r)+a,2Ns+f3{x,da:) 
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(5.13) 
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in which the exphcit form of Z2i<!{s~r)+a,2Ns+(}{^idx) is given by (j3.53p . For example, the first 
term in the last line of (j5.13|) . picking up only the first order in Sk{—dx) and Sk{dx), is given as 

Z2N{s-r)+a,2Ns+a+N{x, Ox) = S2N{s-r)+a+k+l-N{x)S_2Ns-a-k{—x)Sl{ — dx) 

+ S2N{s-r)+a+k-N{x)S_2Ns-~a-~k+l{-x)Si{dx) H 

+ S2N{s-r)+a+k{x)S-2Ns-a-k{—x)SN{ — dx) 

+ S2N{s-r)+a+k-Nix)S-2Ns~a-k+Ni-x)SNidx) -\ • (5-14) 

On the other hand, from (j3.57p . we have first leading terms of the r-function, tm=n{x, g), 
for the LMG model with indices of the Pliicker coordinates, im = A^, + 1,..., = 
— 1, . . . , ii = 1 and etc., in the following form: 



rN{x,g) = v^]'j!f_l'---]l{g) det 
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So{x) Si{x) 















I N,N-1,...,1 ^ \ T , 
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^N,N-i,...,ii9) + VN+i,N::'{l.,ii9) (Siix) + ■■■ + SNix)) + ■■■ 



SN-lix) 
Sn-2{x) 
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(cosf)^ e-i^(^+^) + (cosf) 
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e 2^ 



So(x) 
))- 
(sin I) 



+ 



(5.15) 



The f-dependent HF equation (j3.59p . idvT]\f{x, g{v)} = Hp^-iif{x, dx, g{v)}Tj\f{x, g{v)} , on the 
Tj\[{x, g{v)) may be expected to give a new u-dependent HF solution, where the u-dependent 
HF operator is given by (|5.13p and (j5.14p and the r-function is provided by (j5.15p . 



5.7 Summary and discussions 

In the preceding section, for the LMG model we have an approximate HF operator and a r- 
function up to the first order in Sk{—dx), Sk{dx), Sk{x) and Sk{—x) which should satisfy the 
u-dependent HF equation. Then, we meet inevitably with an interesting and exciting problem 
of solving the u-dependent HF equation on the T^ix, g{v)). After determining HF parameters 9 
and ^jJ self-consistently, a further study should be made to obtain a soliton solution derived 
from a u-dependent Hirota's bilinear equation regarding v as time t and relationship between 
a collective motion and the soliton solution. Such attractive problems have not been treated yet 
and just begin to open. 

Here, we will recur to the representation of the infinite-dimensional LMG model. In the 
previous section, we already have obtained the expression for the LMG Hamiltonian (j5.9p in 
which, however, the commutator term in the first term in the forth line of the equation brings 
an anomaly (an infinitely divergent result) for us, as shown below, 

oo 

y~^(2n -|- 1) ^ oo (anomaly). 

n=0 

For the present, we ought to discard this anomalous term to construct the anomaly-free infinite- 
dimensional LMG model. 



m=0 ^^2m+l 
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Finally, we will point out the possibility of the extension of the present algebra s/(2,C) to 
the af fine Lie algebra A^^^ . Corresponding to an extension of the adopted simple LMG model 
to the so-called coupled LMG model [73\ [71] , we have a very interesting problem of constructing 
an A^^^ LMG model for which the idea given in the paper [51] is considered to be very suggestive 
and useful. According to [51], we can choose the Chevalley basis for A^^^ as follows: 



hi 


= [ei,/i], 


h2 


= [e2,/2], 


ei 






= ^'^2utp*l+2u^ 


fi 




/2 


= ^^'^\+2vi>ly, 



in which the total Hamiltonian of the coupled-LMG system is composed of two LMG model 
Hamitonians }i\ = i/i (J^i^O; -^i;±) and if2 = -f^2(-^2;0) -^2;±) and an interaction term K\-^j^K2 - 
-\-K2-+Ki--.. We will discuss the above A^^'^ LMG model elsewhere using the present infinite- 
dimensional SCF method in r-functional space on Fqo- 



6 Summary and future problems 

In Section [H from algebro-geometric viewpoint, we have given a brief history of microscopic 
understanding of theoretical nuclear physics. It is summarized that we seek for an optimal 
coordinate-system describing dynamics on a group manifold based on a Lie algebra of fermion 
pairs. The TDHF/TDHB are nonlinear dynamics owing to their SCF characters. Seeking for 
collective coordinates in a fully parametrized dynamical system is exactly finding a symmetry of 
an evolution equation in nonlinear dynamics. In differential geometrical approaches for nonlinear 
problems, the integrability conditions are stated as the zero curvature of connection on the 
corresponding Lie groups of systems. Nonlinear evolution equations, e.g., the famous KdV and 
sine/sinh-Gordan equations and etc., come from the well-known Lax equation [70] which arises 
as the zero curvature [29]. These soliton equations describe motions of the tangent space of local 
gauge fields on a time t and a space x, which are Lie group/algebra- valued-equations arising 
from the integrability condition of gauge field with respect to t and x. In the TDHFT/TDHBT, 
the corresponding Lie groups are unitary transformation groups of their ortho-normal bases 
dependent on t but not on x. 

In Section [21 along the Lax form for integrable systems, we have studied essential curvature 
equations to extract collective submanifolds out of the full TDHF/TDHB manifold and shown 
the following: 

(i) Expectation values of the zero curvatures for a state function become a set of equations of 
motion, imposing weak orthogonal conditions among infinitesimal generators, i.e., equations for 
tangent vector fields on the group submanifold. Those of non-zero curvatures become gradients 
of a potential arising from a residual Hamiltonian along collective variables. These quantities 
are expected to give a criterion how the collective submanifold is truncated well. 

(ii) The zero-curvature equation in QPF is nothing but the FRPAEQ imposed by the weak 
orthogonal conditions and has a simple geometrical interpretation: Relative vector fields made 
of the SCF Hamiltonian around each point on an integral curve constitute solutions for the 
FRPA around the same point which is in turn a fixed point in QPF. It means the FRPA is 
a natural extension of the usual RPA for small-amplitude fluctuations around a ground state to 
RPA at any point on the collective submanifold. The enveloping curve, made of a solution of 
the FRPA at each point on an integral curve, becomes another integral curve. The integrability 
condition is the infinitesimal condition to transfer a solution to another solution for the evolution 
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equation. Then the usual RPAEQ becomes nothing but a method of determining an infinitesimal 
transformation of symmetry if fluctuating fields are composed of only normal modes. 

In Section [3l to go beyond a perturbative method with respect to collective variables to ex- 
tract large-amplitude collective motions, we have studied an algebro-geometric relation between 
SCFM and r-FM, method of constructing integrable equations (Hirota's equations) in soliton 
theory. At the beginning, descriptions of dynamical fermion systems in both the methods had 
looked very different manners at first glance. In abstract fermion Fock spaces, each solution 
space of dynamics in both the methods is the corresponding Grassmannian. There is, however, 
a difference between finite-dimensional and infinite-dimensional fermion systems. In spite of 
such a difference, we have aimed at closely connecting the concept of mean-field potential with 
gauge of fermions and at making a role of loop group clear and consequently we have shown the 
relation between both the methods: 

(i) The Pliicker relation on the coset variable becomes analogous to the Hirota's bilinear form. 
The SCFM has been mainly devoted to the construction of boson-coordinate systems rather 
than to the construction of soliton solution by the r-FM. It turns out that both the methods 
are equivalent with each other from the viewpoint of the Pliicker relation or the bilinear identity 
equation defining Grassmannian. 

(ii) The infinite-dimensional fermion operators are introduced through Laurent expansion 
of the finite-dimensional fermions with respect to degrees of freedom of fermions related to 
a V dependent mean-field potential. Inversely, the mean-field potential is attributed to gauges 
of cooperating infinite-dimensional fermions. The construction of fermion operators can be 
contained in that of a Clifford algebra. This fact permits us to introduce an affine KM algebra. 
It means that the usual perturbative method with respect to collective variables with time 
periodicity has implicitly stood on a Groo- Then we rebuilt the u-HFT with the use of the 
affine KM algebra and map it to the corresponding r-functional space. As a result, the u-HFT 
becomes a gauge theory of fermions and the collective motion appears as the motion of fermion 
gauges with a common factor. The physical concept of quasi-particle and vacuum in the SCFM 
on is connected with the Pliicker relations. Extracting sub-group orbits consisting of loop 
paths out of the Grj\f is just the formation of the Hirota's bilinear equation for the reduced KP 
hierarchy to su{N) (c sl{N)). The present theory gives the manifest structure of gauge theory 
of fermions inherent in SCFM and provides a new algebraic tool for microscopic understanding 
of the fermion many-body system. 

(iii) Through the investigation of physical meanings for the infinite-dimensional shift oper- 
ators and the conditions of reduction to sl{N) in r-FM from the loop group viewpoint, it is 
induced that there is the close connection between collective variables and spectral parameter 
in soliton theory and that the algebraic mechanism bringing the physical concept of particle 
and collective motions arises from the reduction from u{N) to su{N) for the u-dependent HF 
Hamiltonian. 

(iv) It must be stressed that though the u-HFT describes a dynamics on real fermion- 
harmonic oscillators, the soliton theory does on complex fermion-harmonic oscillators. This 
suggests us an important task to extend the u-HFT on real space 'su{N) to that on complex 
space sl{N). It gives us a deeper understanding of the concept of quasi-particle energies and 
boson ones, in other words, independent particles and mean-field potential, in a microscopic 
treatment. Recently Wiegmann et al. have developed an approach in which the theory of classi- 
cal integrable systems is applied to studies of ID-fermion systems and the so-called orthogonality 
catastrophe in a Fermi gas. They have introduced a boundary condition changing operator [75j 
but have made no map ctm : : ipiipj Zij{x,dx) (|3.53p contrary to the present u-HFT. 

In Sectional we have given a geometrical aspect of RPAEQ [501 [69] explicit expression 

for the RPAEQ with a normal mode on F^o- We also have argued about the relation between 
a loop collective path and a FRPAEQ. Consequently, the usual perturbative method is shown to 
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be involved in the present method which aims for constructing TDHFT on the affine KM algebra. 
It turns out that the collective submanifold is interpreted as a rotator on a curved surface 
in the Groo- The present theory may lead to multi-circles occurring multiple parameterized 
collective motions. If we could arrive successfully at such a final goal, the present work may give 
us important clues for description of large-amplitude collective motions in nuclei and molecules 
and for construction of multi-dimensional soliton equations |6H [66] since the collective motions 
usually occur in multi-dimensional loop space. 

In Section \5\ as an illustration we have attempted to make a u-HFT approach to an infinite- 
dimensional matrix model extended from the finite-dimensional su{2) LMG model [32]. For 
this aim, we have given an affine KM algebra sl{2,C) (complexification of Hi(2)) to which the 
LMG generators subject and their t representations and the ax mappings for them. Further 
we have introduced infinite-dimensional "particle" and "hole" operators and operators Kq and 
K± defined by the infinite-dimensional "particle-hole" pair operators. Using these operators, 
we have constructed the infinite-dimensional Heisenberg subalgebra of the affine KM algebra 
sl{2,C). Thus the LMG Hamiltonian and its HF Hamiltonian have been represented in terms 
of the Heisenberg basic-elements whose representations are isomorphic to those in the corre- 
sponding boson space. They have been expressed in terms of infinite numbers of the variables 
Xfc and the derivatives dx,. through the Schur polynomials Sk{x). Further we have obtained an 
approximate HF operator and a r-function up to the first order in Sk{—dx), Sk{dx), Sk{x) and 
Sk{-x). 

In Appendices, we have given the infinite-dimensional representation of SU{2N)oo transfor- 
mation of the "particle" and "hole" operators. The expression for t rep of (Y„(2i+i) + Y(2i+i)), 
i.e., <7y„{2j+i)+Y(2,+i) {^) lias been first given in terms of the Bessel functions. We have also shown 
an explicit expression for Pliiker coordinate and calculated a quantity, det(l7v +P^p), in terms 
of the Schur polynomials. 

Finally intimate relation of SCFT to soliton theory has been shown to come from ways of 
constructing a closed system of solution spaces. The ordinary SCFM has been almost devoted 
to approach cooperative phenomena in finite fermion systems. We must contrive construction 
of the optimal coordinate-system on the group manifold. For this purpose the relation between 
the boson expansion method for finite fermion systems and the r-FM for infinite ones should 
be intensively investigated to clarify algebro-geometric structures of integrable systems. Such 
algebro-geometric approach will make a bridge between finite fermion systems and infinite ones. 
Various physical concepts and mathematical methods will work well also in the infinite ones. 
The SCFM based on global symmetry should be much improved noticing local symmetry of the 
infinite ones and then may open a new area in vigorous pursuit of wider fields of physics. 

We have many future problems in connection with the above discussions, which are itemized 
as follows: 

(i) To study the relation between the quantity of non-zero curvature and the collectivity: It is 
interesting to study the relation using the simple LMG model, which leads to an investigation of 
the effective condition for the collective submanifold extracted by the zero-curvature equation. 
Temporarily digressing from the integrability condition, adopting the Bethe anzatz (BA) we 
have obtained exact solutions for the LMG model solving the BA equation [76]. Contrary 
to Pan and Draayer's work [63] and our previous works [77], we do not use any bosonization 
nor infinite-dimensional techniques and hence have no restrictions on interaction-strengths of 
LMG Hamiltonian. Considering the advantage of the integrability condition, the famous Gaudin 
model plays an important role to solve effectively the BA equation [78]. From the loop group 
viewpoint, as shown by Sklyanin, with the use of the exactly-solvable Gaudin model obeying the 
Gaudin algebra, an exponential generating function of correlators is obtained from the Gauss 
decomposition for sl{2, C) loop algebra, which gives correlators including the Richardson-Gaudin 
determinant formula for the Bethe eigen- function [79] . A generalization of the Gaudin algebra is 
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given by Ortiz et al. [80]. These works may have an intimate relation with our recent work [35] 
and the present work. 

(ii) To clarify the explicit relation between spectral parameter and collective variable and 
the physical concept of the geometrical connection: The spectral parameter of the iso-spectral 
equation in soliton theory and the collective variable in SCFM, though showing different aspects 
at a glance, work as scaling parameters on S^. The former relates to a scaling by analytical 
continuation of S^, i.e., z. The latter makes roles of deformation parameters of loop paths 
in GrM- 

(iii) To study the relation between weak boson operators and boson mapping operators, i.e., 
the shift operators in r-FM: The generators for collective variables in Fqo can never take exact 
boson commutation relations because of the finite-dimensional matrices. 

(iv) To study a relation hidden behind gauge of state functions and construction of fermion 
pairs: In the usual algebraic treatment of fermion many-body systems, we assign an abstract 
number to each of the set of quantal numbers and let their fermions make the Lie algebras 
{u{N), so{2N), so{2N + 1) and etc.). For the pair-constructions we have an interpreta- 
tion as classifications of Laurent spectra in the infinite-dimensional fermions, although we did 
not manifestly state. On the other hand, as well known, electron spin can be described as 
a geometrical phase of gauge with the help of Mobius band. Then we inversely start from 
fermions with the abstract numbers and through any way we could return to the original 
fermions with the physical quantal numbers. We think it not so wrong to attempt to under- 
stand quantal numbers as the geometrical attribute of the Grassmannian made of the abstract 
fermions. 

(v) To establish mathematical tools to obtain subgroup-orbits of loop paths in Gtm, basing 
on the Pliicker relations on S^, i.e., soliton equations. This problem is most fundamental to 
solve the new theory: For this aim we must know sub-group orbits or corresponding sub-Lie 
algebras and establish mathematical tools to extract them out of Grjv/. In this concern, we are 
intensively interested in the algebraic mechanism for spontaneous decision of a fixed point and 
a collective submanifold around the point. 

(vi) To study a relation between nonlinear superposition principle in soliton theory and 
generator coordinate method (GCM) in SCFM: The GCM may provide a superposition principle 
on a nonlinear space [811 116j . Standing on the viewpoint of local symmetry of infinite fermion 
systems behind global symmetry of finite ones, we might reconstructed the GCM and nonlinear 
superposition methods using the infinite-dimensional shift operators. What relation does exist 
between the construction of exact solutions based on the idea of the imbricate series in soliton 
equations ^2] and resonating mean-field theories [83j ? 

(vii) To study why soliton solutions for classical wave equations show fermion-like behaviours 
in quantum dynamics and about what symmetries are hidden in soliton equations. To both 
the questions suggested by Tajiri et al., we cannot give a satisfactory answer yet within the 
present framework. Because both the methods are a priori based on the fermion system from 
the outset. That is to say: SCFM describes a quasi classical dynamics on Grassmannian (S- 
det orbit) which is induced owing to the anti-commutative property of fermions. On the other 
hand, r-FM also uses the fermions to explain Grassmannian of solution space which reflects the 
fermion-like behaviours of soliton solutions. Therefore we should study further why extraction 
of soliton equations out of classical wave equations brings out Grassmannian. In the reductive 
perturbation method [H^, soliton equations appear as the symmetry space in a classical wave 
equation with respect to a transformation and a scaling transformation with a common parame- 
ter for independent and dependent variables. We can see a relation between the power exponents 
of the parameter and the degree of shift operators. At any point on extracting way of soliton 
equations from classical wave equations did we introduce the "anti-commutative character", in 
other words, did we introduce the structure of Grassmannian? 
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A Coset variables 

Following [TB], let us introduce the triangular matrix functions S{(), C{C) and C{C,) defined as 



ACQ 



s(c) = (s„(c)) = g(-i) V+ i)r 

OO 

OO ^ 



k=l 



which have the properties analogous to the usual triangular functions 

C\C) + SHC)S{C) = Im, C\C) + SiOS^iC) = In^m, S{C)CiC) = (7(05(0. 

Then the matrix p in (j3.3p is defined as p = (pia) = S{C)C~^{C)- The matrix g in p.3p is 
decomposed as g = gc^Qw using the matrices given below 



5C 



C{C) -5(0t 
s{0 c{0 



1 


" 


-cn 


7 = 


. C 






w = e' 



w = e' 



7 



T] 
7? 



V 



-V, 



9w = e 



,7 



w 
w 



-V, 



where C is a (A^ - M) x M matrix (Qa) and 7] and f] are M x M and (TV - M) x (A^ - M) 
anti-hermitian matrices (r/afe) and (??ij), respectively. The indices i and a denote unoccupied 
(M + 1 ~ iV) states and occupied (1~ M) states, respectively. 

B Properties of the differential operator eia acting on $m,m 

Let us introduce a differential operator and a function corresponding to the so-called particle- 
hole operator and the free particle vacuum in the physical fermion space, respectively as 



PibPja 



d 



+ 



d 



d 



dpjb dp*^ 2 dr 



$M,A./(p,/,r) = [det{l+p^p)]-2e'\ (B.l) 

By using the formula for the differential of a determinant, we can easily calculate as 
d 



dpjb 
d 



[det(l +ptp)]-^ = -ip*J(l +ptp)-i]T[det(l 
[det(l +ptp)]-| = -Ipiaiil +ptp)-i]^^[det(l 



(B.2) 



^^^^ 

where denotes an inverse matrix of a. Then from equations (jB.lh and (IB.2P we get 

eta^MM = \j2PibPja ■ Pjc[(l +p'^P)~^]lb + +pM~"^]da + \Pia^ '^M,M 

= i {p{l+p^py^{l+p^p)}ia + \pia ^M,M =Pia^MM^ 
4aPjb = -PibPja + Pjbe*a- 

Thus we can prove 

eta^M,M = Pia^M,M^ [f^ia^Pjb] = -PibPja- 

These are just the relation given in the first line of (j3.11|) . The other relations in (|3.11|) can be 
proved in a similar manner |37j . 
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C Affine Kac— Moody algebra 



According to Kac and Raina [12] , let gl{N) be the Lie algebra of all x matrices with complex 
entries acting in C'^ and let C[2;, z"^] be the ring of Laurent polynomials in indeterminate z 
and The loop algebra [36] gl{N){Du{N)) is defined as gl{N){C[z, z^^]) (Dn(A)(C[z, z"!])), 
i.e., as the complex Lie algebra of A x A matrices with Laurent polynomials as entries. An 
element of gl{N) is given in the form 



r 

Z ttr 



{ar £ gl{N)). 



(C.l) 



As was pointed out by Kac and Raina, we regard the fermion pair- and single-operators as 

^ eai3 (1 < a, /3 < A), 
Ua (1 < a < A), Co 



/^t I — V 



(1 < a < A), 



where the matrix elements ea/s are equal to 1 in the (a, /5) entry and elsewhere. The matrix 
elements e^^ form a basis of gl{N). The components of A x 1 column vectors Ua are equal to 1 
in the a-th row and elsewhere. They span the vector space in which the gl{N) acts. The 

symbol T means the transpose of a vector or a matrix. The matrices eap{r) = z^Cap constitute 
a basis of gl{N). The loop algebra gl{N) acts in the vector space C[z, z~^]'^ consisting of A x 1 
column vectors with the Laurent polynomials in z and z~^ as entries. The Lie bracket on gl{N) 
is the commutator 



[e„/3(r), e^s{s)] = Sp^Casir + s) - dase^pir + s) 
The vectors defined as 



(C.2) 



^Nr+a 



—r 

Z Ur 



T r 

U^Z , 



also form a basis of the vector space C[z, z ] indexed by Z and its dual space. The {fj | i = 
Nr + a € Z} is given by the column vector with 1 as the i-th row and elsewhere. Thus it is 
possible to identify C[z, z"^]^ with C°°. The relation ea(}{r)vNii^ = VN{s-r)+a is easily derived. 

For a{z) S gl{N) we denote the corresponding matrix in aoo by r{a{z)}. Then we deduce 
a matrix representation for T:{eaf3{r)} in cloo as 



'^{eQ/3(r)} = j^Af(s-r.)+a,Afs+/3- 



Eij G 
algebra Hoo 



have 1 as the entry and elsewhere and form a glc 



for \i - j\ > N}. 



(C.3) 

Suppose a bigger 
(C.4) 



There exists such an N satisfying the above condition. The corresponding matrix of the a{z) 
in (jC.ip in has an infinite A periodic sequence of block form 



T{a(z)} 



a_i oo ai 

a_i oq ai 

a_i Oo oi 



(C.5) 
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We regard (jC.SP as a representation of the matrix Cqo in which elements on each diagonal parallel 
to the principal diagonal form a periodic sequence with period N. Let X{k) = z^X {X G gl{N)) 
be an element of gl{N). Define an antilinear anti-involution to on gl{N) by Lo[X{k)] = z~^X'^ . 
Then in the a^o we get x{oj[X{k)]] = Tt{X(A:)}. 

Using the fundamental idea of the Dirac theory [40], we define the vacuum state in which 
the state labeled by the "Laurent spectrum" with positive energy is empty but all the negative 
energy states labeled by the "Laurent spectra" are occupied. Denoting an exterior product of 
vectors as A, a perfect vacuum ^'o and a reference vacuum ^! m are expressed as 

$0 = l/Q A Z^_l A Z^_2 A • • • , ^A/ = l^M A VM~l A A • • • . 

Let the space V = (BiezCvi and its dual V* = ©jg^Ci^'* be an infinite-dimensional complex 
space with a basis {vi, v* | i, j G giving a linear functional v* on V by = 5ij G Z). 

Each 1/ £ V and v* G V* defines a wedging operator v and a contracting operator u* on the -Fqo 
as 

v{vi^ A f j2 A • • • ) = ^ A z^i-^ A z^i2 A • • • , 

v*{vi^ A A • • • ) = ^*{vii)^i2 A 1^13 A A v*{vi^)vi^ A A f A • • • 

+ V* {vi.j)vi-^ ^ ^ ^ ■ ■ ■ . (C.6) 

Then the operators {vi, u* \i,j G Z} generate a Clifford algebra 

K,i;,} = 6ij, {ulu*} = = 0. (C.7) 

Thus the anti-commutation relations lead us to the identification of the new fermion annhilation- 
creation operators (I3.23P with the present operators at a pointwise z with |z| = 1 as 

Using the above identification, the corresponding perfect vacuum and reference vacuum, it can 
be shown that 

^-o^lVac), I'm ^ |M) = VM---^/'i|Vac), 

ipNr+a\ya'C) = 0, (Vac|-0?^^+Q, = (r < -1), 
V'^^+^lVac) = 0, (Vac|^7vr+a = (r>0). 

From (|C.3p . (|C.6p . (jC.Sp and the Chfford algebra (|C.7p we obtain representations 

T^{ea/3(r)} = ^ S^(s_r)+a,Ar^ - ^ ^'Ar{s-r)+ai^iVs+/3 - ^Af(s-r)+aV'Ars+/3- (C.9) 

sgz sez sez 

Construction of the SCFM on the infinite-dimensional Fock space is an interesting illustrative 
problem in order to clarify an algebraic structure among the original fermionic field, the vacuum 
field defined in the SCF potential and the bosonic field associated with Laurent spectra. 

It is well known that the operator of ()C.9p acting in the F^o has in general an "anomaly". 
To avoid the "anomaly" , we had better use either of the normal-ordered products given below 

: ENr+a,Ns+l3 '■= Ei^r+a,Ns+P — ^a^^rs {s < 0), 

: UNr+a^^Ns+fd ■~ ^Nr+ai^Ns+fi ~ ^ajS^rs {s < 0), 

: Ip^r+ai^hs+P '■= '^pNr+a'fphs+P " ^afi^rs (s < 0). (C.IO) 
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We define important shift operators A-,- and A_j (j G Z_|_) in the sohton theory on a group as 



1 



z 



1 




A_,- = T 



•• 

1 ■■ 







1 



(C.ll) 



which are images under t as 



A 



■0 



a=l 



'N-l 



clca+l+clci^ 



a=l 



Through an easy calculation with the use of ()C.9P and (jC.lOp . the 2-cocyde a on a^o ()C.4p 
induces the following 2-cocycle on a pair of basis elements of gl^: 



a{T:{eai3{k)} ,r{e^s{l)}) = ^aS^iSj^k+ifi ■ k, 
where the Kac-Peterson 2-cocycle a [85l[l8l[5l] is defined as 



(C.12) 



a{Eij, Eki) = 6jkdii (i < 0) - 6ii6kj {j < 0) 



1, for j = k, j > 1, i = I, i < 0, 
-1, for i = /, i > 1, j = k, j < 0, 
0, otherwise. 



Then we have for the shift operators (|C.lip 

a(Afc, A/) = 5k+ifl ■ k. 
For any elements a{z) and b{z) in the gl{N) the formula (|C.12p is written as 

a{x{a{z)],r{h{z)]) = ResoTr a'(2)6(z), 

where a' {z) is the derivative of a with respect to z and Reso is the residue at z = 0, i.e., the 
coefficient of z~^. Investigation of the highest weight representation oi gl{N) leads to its central 
extension gl{N) = glj^ + C • c in which general elements a{z) and h{z) and center C • c satisfy 
the KM brackets 



[a{z),c]KM = 0, 

[a{z),h{z)]]^M = [a{z),h{z)] + {ResoTr a'(z)6(2;)} • c. 



(C.13) 



The Lie algebra gl{N) is called the affine KM algebra associated with the Lie algebra gl{N). 
For simplicity consider the level one case, c\M) = 1 • \M). Using the one-level formula it is 
possible to rewrite (|C.13p as 



Xa=Xa + C-C, 
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[Xa, c]kM = 0, [Xa, Xf 



b\KM 



X[a,b] + a(«, b) • c. 



The matrices a and 6 are any elements of (IC.Sp . The [a,b] denotes the Lie bracket of their 
matrices. The T-rep of a{z) is given as 

T{a(z)} = ^ y^^{ar)N{s-r)+a,Ns+l3''pN(s-r)+a'4'Ns+g^ 

rez sez 

from which the commutator and the 2-cocycle a between r{a{z)} and 'x{b{z)} are calculated as 

Ma{z)},r{b{z)}] = ^ ^{[ar,as])N(t-r-s)+a,Ns+gi^N(s-r)+a1p*Nt+p^ 

a{r{a{z)},r{b{z)}) = ^ r TV(a,6_^). 



D Schur polynomials and r-function 

The Schur polynomials Sk{x) belonging to C{xi,X2, • • • ) are defined by the generating function 



exp ^ Xkp'' = ^ Sk{x)p^ 



(D.l) 



k>0 



The Schur polynomials and their associated recursion formulas play important roles to evaluate 
matrix elements between two number-projected Hartree-Bogoliubov states [86]. The Schur 
polynomial is related to a symmetric function /i^. ^ hf^p^ = n^(l — eip)~^. Then the Schur 



A:>0 



polynomial Sk{x) is written as 

Skix) = /ifc(ei,e2, . . . ,eAr), 
and is given explicitly as 

S'o(x) = 1, Si{x) = xi, 



1 J 



ei + €■!, + ■■■ + e 



N 



S2{x) = X2 + \xl, 



S^{x) = X3 + XiX2 + lx\, 



S'4(x) = X4 + X1X3 + \X2 + \x\x2 + ^x\, 



We here construct an explicit form of the polynomials in the bosonic Fock space which under 
the (TM corresponds to the finite monomials of the F^^^\ It is given by 



where iM > Uf-i > • • • > ^i- The Schur polynomial Sx{x) = 5'Ai,A2,A3,...,Afc(2;) is given as 



(D.2) 



Sx,,X2M,...,Xf,{x) = det 



Sxiix) Sx-,+iix) Sx^+2ix) 

Sx2-i{x) Sx2{x) Sx2+i{x) 
Sxs-2{x) Sx^^iix) Sx-^{x) 

Sxk+i~k{x) Sx^+2~k{x) Sx^+3-k{x) 



Sxi+k^iix) 

Sx2+k-2{x) 
Sxs+k-six) 



where A denotes a partition Par A = {Ai > A2 > • • • > Afc > 0} [481 149j . As for the contravariant 
hermitian form (j3.52p in they form an orthogonal basis {Sx\Sf^) = Sxfi- For the special 
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partition Par A = 1^^ = {Ai = 1,A2 = IjAs = 1,A4 = 1,...,Am = 1}, i-e., the completely 
anti-symmetric Young diagram, 



SiM (x) = det 



Si{x) S2{x) Ssix) S^ix) 
1 Si{x) S2{x) Ss{x) 
1 Si{x) 52 (x) 



















Sm{x) 
Sm-i{x) 
Sm-2{x) 

Si{x) 



{-lf'SM{-x). 



A group orbit of the highest weight vector \M) under the action U{g) is mapped to a space 
of r-function TM{x,g) = {M\e^'^^'^U{g)\M) . Using (|D.2|1 . the Schur-polynomial expression 
for T-function is given by (j3.56p . Noting the relation {S\\S^) = Sx{dx)Sfj_{x)\x=Q = 5\^, 



M,M-1,---,1 



'^ "\^{g) is obtained from TM{x,g) as 



M,M-1, dots,!. 



"*M,»A/-i,.''°il'^(5') = '^iM'^M-i',--'!n"L5'(a^)}U=0 = Sij^.j_M,iM-i-{M-l),...,h{dx)TM{x, g)\x=o, 

oo 

v{g{x)} ^ e^(^)i;(g)e-^(^), g,j{x) = Y,Sk-^{x)gkJ{x). 

k>i 

The coefficient vf^'^~^'"''\ (q) is called the Pliicker coordinate 1871 . In the soliton theory on 

a group the PI 
KP hierarchy. 



a group the Pliicker relations among vf^j'fj^j^l '\-i^{9{x)} correspond to the Hirota's forms in the 



E Hirota's bilinear equation 

According to Kac [18], a vector |tj\/) G belongs to the group orbit of the highest weight 

vector (the vacuum) if and only if it satisfies the bilinear identity equation 

Y,i^i\TM)^ri\rM) = ^ \TM) = U{g)\M), g e GL{oo). 



This identity is cast into an infinite set of nonlinear differential equations for TM(x)={M\e^^^^\TM] 

= -L / ^(M + l\e"(-')^^,(p)\TM) ®{M- l|e^("")^*(p)|TM) 
zvrz J p 




It is described by the Hirota's bilinear differential operator as 



P{D)f -g^P (^^, • • • J (fix + y) ■ g{x - y))\y=^, 

where P{D) is a polynomial in D = {Di,D2, . . . ). Note that Pf • / = if and only if P{D) = 
—P{—D). Defining new variables x and y by the relations x' = x — y and x" = x + y, with the 
help of (jP.ip and the notation D = {Di, ^D2, . . .), (jE.ip is brought to the form 

S,i-2y)Sj+i{D) exp I y^D, | tm(x) • tm{x) = 0. (E.2) 

j>0 \s>l J 

If we expand ()E.2p into a multiple Taylor series of variables yi,y2, ■ ■ ■ and make each coefficient 
of this series vanishing, we get an infinite set of nonlinear partial differential equation for the 
KP hierarchy. 



d d 
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F Calculation of commutators and 2-cocycles among 
operators K 

Commutators and 2-cocycles among operators i^o(^) and K±{r) are calculated using (|C.2p 
and ()C.12p as follows: 



[Ko{r),K±{s)] 



\ Jap V j2N 



'P \ f2N 



{6i3^ea5{r + s) - Sase^pir + s)} 



e2N 



T2N 



eas{r + s) 



aS 



hN J 



\h2N, 



e2N 



T2N 



dapir + s) 



aP 



e2N 

-f2N 



eap{r + s) = ±K±{r + s), 



aP 



a{Ko{r),K±{s)} = a{x{Ko{r)},x{K±{s)}} 

fh2N) Meap{r)}, 

/ at) 



\h2N 



' aP 

e2N 
aP V f2N 



f2N 



r{e^s{s)} 



7<5 



■yS 



iXr ( h2N 



e2N 



2N 



Sr+s,0 • r = 0, 



[i^+(r), A'_(s)] = {e2N)apeap{r),{f2N)-y5e^5{s) 

= {e2N)aP X {f2N)yS p-f^asir + s) - SaS&ypir + s)} 
= {e2Nf2N)aSea5{r + s) - {f2Ne2N)ypeyp{r + s) 

[e2N,f2N]) eQ/3(r + s) = 2 (i/i2Ar) eap{r + s) = 2Ko{r + s 



aP 



aP 



a{K^{r),K.{s)} = a{r{K+{r)},r{K.{s)}} 

= a^{e2N)apT:{eapir)}, {f2N)'yS'^{e-yS{s)}^ 

= {e2N)ap{f2N)-ySSaSSpjSr+s,0 ' r = {e2N f2N) ^r+s,0 ' r = Nr6r+s,0- 

Commutators and 2-cocycles among operators Xk and Yk (fc = 2r + 1) are calculated as follows: 

[Xk,Xi] = [Yk,Yi] = 0, 



a 



XA /X, 



Y, r\ Y 



a 



a <r 



Xi. \ / X i 



Yi 



{k = 2r + 1, / = 2s + 1) 

|(e2Ar)«/3T{eo/3(r)} ± if2N)apT:{eapir + 1)} , 
(e2Ar)75T{e^5(s)} ± {f2N)yS'^ {e^sis + 

{e2N)ap{e2N)-y5Sa5Sp'ySr+s,0 ' f 

± {e2N)ap{f2N)'y5SaSSpjSr+s+l,0 ' r ± {f2N)apie2N)'ySSaSSpjSr+l+s,0 ' + 1) 
+ if2N)ap{f2N)'ySSa5Sp'ySr+l+s+l,0 • (r + 1) 
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[-'^fc,^/] = [{e2N)al3ea(3{r) + U2N)al3eaf}{r + 1), (e2Ar )'y5e^5 (s) - {f2N)^Se.^s{s + l)\ 

= [{e.2N)af3 X {e2N)-i5\ {^p^^easir + s) - daseji3{r + s)} 

- [ie2N)al3 X {f2N)-f5] {Si3^eas{r + s + 1) - 6aSe-yi3{r + s + 1)} 
+ [{f2N)af3 X {e2N)-y5] {^f^^easir + 1 + s) - dase-yfsir + 1 + s)} 

- [(/2Af)a/3 X {f2N)^5] {^f^^easir + 1 + s + 1) - 6ase^p{r + 1 + s + 1)} 

= {e2Ne2N)a5eaS{r + s) - {e2Ne2N)-yl3e-yl3{r + s) 

- {e2Nf2N)a5ea5{r + S + 1) + (/27Ve27v)^/3e^_a(r + s + 1) 
+ {f2Ne2N)a5ea5ir + I + s) - ie2N f2N)'yf3e^p{r + 1 + s) 

- if2Nf2N)a5ea5{r + 1 + S + 1) + {f2N f2N)^pe^fi{'r + 1 + S + 1) 
= -([e27V,/2Ar])„^e„/3(r + S + 1) + ([/2Af,e2Ar])^^e„/3(r + s + 1) 

= 2{-h2N)apeaf3{r + S + 1) = 2y2r+l+2s+l = 21^+;, 

f3{r) + (/2Af) MO) 

= [(e27v)a/3, (-^2Af)75] {'5/37ea5('^ + s) " ^aSe^/sir + s)} 

+ [(/2Af)a/3, i-h2N)-y5] {S/S^easir + I + s) - SaSe^/sir + 1 + s)} 

= — (e2Af/l2Af) (r + s) + (/l2Afe2Ar)7/3e^/3(r + s) 

- {f2Nh2N)a&ea5{r + S + I) + {h2N f2N)^l3e-yf3ir + s + I) 

= {[h2N,e])^peap{r + s) + ( [/i, /2Ar] )^^e„;3(r + s + 1) 

= 2{(e2Ar)a/3eQ,/3(r + s) - {f2N)af3eal3{r + S + 1)} = 2y2{r+s)+l = 21fc+25, 

a(Xfe,y2s) = a{T(Xfe),T(y2s)} 

= a{{e2N)ai3T:{ea/3ir)} + (/2Ar)a/3T{eQ/3(r + 1)} , i-h2N)'y5T:{e^5{s)} } 

= {e2N)al3{ — h2N)^5Sa5Sl3'r^r+s,0 ■ f 

+ {f2N)ap{-h2N)^55a5^f3'y^r+l+sfi • (r + 1) = 0, 
[Yfc, 11] = [(e27v)a/3ea/3(?^) " {f2N)al3eaf3{r + 1), (-/l27v)7<5e75(s)] (/c = 2r + 1, / = 2s) 
= [(e2Af)a/3 X (-/l27v)75] {(5/37ea5(r + s) - (5a5e^/3(r + s)} 

- [(/27v)a/3 X {-h2N)^5] {SfS^Casir + I + s) - 5a&e^p{r + 1 + s)} 

= ([/i27V,e2Ar])„^e„/3(r + s) - ([/i2Ar,/2Ar])^^eo/3(r + s + 1) 

= 2{(e2Ar)„/3ea/3(r + s) + {f2N)apea(3{r + s + 1)} = 2Xfc+i, 

1 



^ Y2r.-^Y2s 



N 



N 



2r5, 



1 



d 



2r 



'N 

a(Xfc,y2s) =0, 



dy2r 



1 



a 



N 





1 








l2. 






KM 



Yn 



2r-> , , , 

^ 2ry2r, 
2 1 



2s 



2r ■ 6r+sfi ■ c, 



KM 



N VN 



Yk+2s • Cj 



(k^-l), 



1 



Yk, 



1 



Iz 



1 



KM 



NVN 



Xk+i ■ c. 
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G Sum-rules for 2{Yi + Y_i) 

First we give the vertex operator for T{p) in terms of the Schur polynomials 
r(p) = expi Yl P™(2x™)>exp| ^ V (-29„) 

oo 
n=0 

According to Kac [58], we have 

5^p-'y, - i(r(p) - 1), 

oo 

2Yi + 6io = Y, Sn-i{2x)Sn{-2d) {xn = dn = for even n). (G.l) 

n=0 

Then using (1G.1|) . the sum-rules for 2{Yi + y„j) are derived as follows: 

2Yo = Si{2x)Si{-2d) + S2{2x)S2{-2d) + 53(2x)53(-2a) 

+ 54(22;)54(-2a) + 55(2x)55(-2a) + Se{2x)S6{-2d) + ■■■ , 

oo 

2(y_2 +Y2) = J2 iSn+2{2x) + Sn-2{2x)) Sn{-2d) 
n=0 

= S2{2x) + 52(-29) + 53(2x)5i(-2a) + 5i(2x)53(-2a) + Si{2x)S2{-2d) 
+ 52(2x)S4(-2a) + S5(2x)53(-2a) + Sii{2x)Si{-2d) + • • • , 

oo 

2(y_4 + Yi) = Y, (Sn+4(2x) + S„_4(2x)) 5„(-29) 

n=0 

= S4(2x) + 54(-2a) + S^{2x)Si{-2d) + 56(2x)52(-2a) 
+ S^{2x)S:i{-2^) + S8(2x)54(-2a) + • • • , 

oo 

2(y_i + Yi) = (S„+i(2x) + S„_i(2x)) Sn{-2d) 

n=0 

= Si(2x) + Si{-2d) + 52(22;)5i(-2a) + 5i(2x)52(-25) 
+ 53(2x)S2(-2a) + S2(2x)53(-2a) + 54(2x)S3(-2a) + • • • , 

oo 

2(y_3 + Y^) = Y, (Sn+3(2x) + S„_3(2x)) Sn{-2d) 

n=0 

= S3(2x) + 53(-29) + 54(2x)5i(-2a) + 5i(2x)54(-29) 
+ 55(2x)S2(-2a) + S6(2x)53(-2a) + 57(2x)S4(-2a) + • • • , 

oo 

2(y_5 +Y^)=Y, (Sn+5(2x) + S„_5(2x)) Sn{-2d) 
n=0 

= S5(2x) + SQ{2x)Si{-2d) + S7(2x)52(-2a) 
+ 58(2x)S3(-2a) + S9(2x)54(-2a) + • • • , 
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H Expression for fl'Y_(2i+i)+Y(2i+i)(^) terms of Bessel functions 

Using the formula for (y_(2i+i) + Y,^2i+i)), 

y_(2j+l) + Y{2i+1) = (e2Ar)a/3eo/3(-i - 1) - {f2N)cxpeal3{-i) 
+ (e2Ar)a/3e«/3(i) - if2N)af3eal3ii + 1), 



e2NZ 



-i I 



f2NZ ' + e2Nz' - f2Nz'^= e2N{z ' + z') - f2N{z ' + Z 



Z^Hz^'e2N-f2N){l + z''^') 



2i+ll 



(H.l) 



we express straightforwardly the flY_(2,+i)+Y{2i+i) (-2), T rep of (jH.ip . in terms of the Bessel 
functions as, 



{2i+l)-|-I{2i+l) 



+Y(2.+ 1) (^) = T 





'In 


■ 









Itv. 



^-(2m)(l + ^2*+l)2.i^ 

(z-(2m))2(i+^2m)4.i^ 

(z-(2*+l))2(l + z2m)4.i^ 



+ 



^2i+U4 ^ 



+ (1)^ z-(.-^e2iV - f2N){l + z'^^') + (f )^ (-,-3-i)(,-ie2^ _ + ,2m)3 

+ir ^^^'-'(^-^62^ - /2iv)(l + z'^-^'f + ■■■] 

{(f )' Z"(2m)(i + ^2.+!) + ^ (1)3 (_^_(2^+l))2(^ ^ ^2^+1)3 



I -^27V 



} f2N 



a(3 



= J^(-1)V2. (f) (i27v)a/3T{e,^((2i+l)r) + e,/3(-(2i + l)r)} - Jq (f) (l27v)a/3T {e„^(0)} 
r>0 

+ ^(-l)V2r+i (f ) (e27v)a/3T{e„;3((2i + l)r + i) + e„^(-(2i + l)r - {i + 1))} 

r>0 

- Y,{-iyj2r+l (f ) (/2iv)a/3T{e«;3((2i + l)r + i + 1) + e^fi{-{2i + l)r - i)} , 



r>0 

where we have used the relations 

5:(-l)V2.(f){.(2^+^>+, 



-(2i+l)r 



} = E(-i)^^2. (I) 



r>0 



5^(-l)V2.+l (f ) {^(2^^+^)^+^ + = ^(-1)V2.+1 



^(2i+l)r+i 



r>0 



^(-1)V2.+1 (f ) {.(2^+1)^+^+1 + .-(2^+1)-^} = ^(_l)V2.+l (f ) Z 



{2i+l)r~i 



r>0 



Finally we can express the Qy 



5Y_(2,+i)+Y(2,+i)(2) = t 



(2i + l)+Y{2i+l) 



z as 



^(_1)V2. (f ) z(2^+l>/2^ 
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+ J](-1)V2.+1 



e\ J2i+ 



r+1 ( 2 j ^ 



^ ' J2N 



^(-1)V2. (f ) (/27v)a/3T{e„/3((2z + l)r)} 



rG2 



+ ^(-l)V2,+i (f ) [{e2N)apr{eap{{2i + l)r + i)} - {f2N)ai3r{eai3i-{2i + l)r - i)} 
This is the first expression which we give in terms of the Bessel functions. 



I Properties of SU{2N)oc transformation matrix 

In (|5.10p the 2A^-dimensional matrix g^^ is expressed as 



hN+i[ -yj h2N + ^ 



2 



h2N + 



3! 



^2N + 



1 



''2N 



1 



2 1 



1 / , 1 / 9J, 



i27V 



cos ( — l2Ar + « sm — h2N = „ -, 

2 / \ I J 62. 1^ 



(I.l) 



Adjoint actions K±{6) for ipNr+a are computed as 

[K+{9),tpNr+a\ = ^ V'Ar(s-r)+/3 {-dre2N)pa ' 

[^_(6l),V'Arr+a] = J]] V'iV(s-r)+/3 ( " ^-r Aiv) , (1-2) 

from which we obtain 

[K+{0) - K_{e),i;Nr+a] = ^ ^Ar(s-r)+/3 {-Or) {e2N - ^n) {0-r = Or)- 

We also have the relations 

{e2N — /2Ar)^ = —I2N, {e2N — /27v)'^ = —{e2N — f2N), 

{e2N — /2Ar)^ = hN, (1-3) 

Using p.2p and (jl.sp . the transformed il^Nr+aide) is shown to be (jS.lip . The 2A^-dimensional 
matrix gg^ is expressed as 

ffe, = i27V + (^-y^ (e27V - Atv) + ^ (^-y^ (e27V - /2Jv)^ + ^ (^-y^ (e2JV - Ajv)^ 



l-^(-|)^ + ^(-|)^ + ---K2^ 



(1.4) 
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+ 



cos I -— I • l2Ar + sin ( -y ) • {e2N - f2N) 



cos 



|)-l;v -sin(|)-ljv 
sin(%)-i^ cos(%)-iiv 



Using (jni) and (Ol) . we obtain an explicit expression for tlie block matrix Qu^t^s of (2A^)c 
transformation-matrix appeared in ()5.12l) as 



-t — 



). (1.5) 



cos (I) e-^('/'"+^=) . 1^ - sin (|) e-sCV---^^) . 
sin (I) e5(^"-^-) • Ijv cos (|) e5(^-+^^) • 1^ 

The transformation (j5.12p is rewritten under the change of index of the block matrix as follows: 



i^Nr+aig) = X] I '4'N{r-s-t-u)+f3 ) {gu,t,s)l3a {gu-t,s = gu,t,s) 



— 2^VN{r-s)+l3 I 2^ 9u,t,s-t-u 



13a 



t,s—t~u 



Noting the definition {g)Nr+a,Ns+f3 = {gs-r)ai3 in the second equation of (j3.28p and again chan- 
ging the index of the block matrix, the above transformation is further rewritten as 

1pNr+a{g) = ■ll^Ns+p{g)Ns+f5,Nr+a = Y^ IpNs+fsjgr-s) Pa = Y, '>pN(r-s)+pigs) f3a 
sel. seZ sez 



'4'N{r-s)+l3 j X] gu,t,s-u-t 



t^t~u 



Pa 



'4>N{r-s)+P j X] 9u,t-u. 



s-t 



Pa 



Then the matrix element of g is also represented as 



(g) Nr-\-a,N s-\-P I / ^ gu^t^s—t—u—r 
, t,u£'Z 



^ ^ gu,t—u. 



s—t—r 



aP 



aP 



J Explicit expression for Pliicker coordinate and calculation of 
det(l7v -\-p^p) in terms of Schur polynomials for LMG model 



From ()3.8p , the coset variable p is expressed in terms of Pliicker coordinates as 

l,...,a,...,M / 

P={Pia) = {[S{C)C~\C)]i 



^i;.'.'.'2'.'.;M (5c) 



1 1,...,M/ N I 

Using explicit expression for SU{2N) matrix of the Thouless transformation (j5.2p . the matrix 
representations of the Pliicker coordinates '^^(^f ) and "'^/(^c) ^^'^ given simply as 



l,...,a,...,M I 



^i,'...,7,'...,M (fc) = sin I • Iat, ^'l,'.'.'.',Af (5c) = cos I • In (M = iV). 
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Then we have p = tan | • Iat, which is identical with the expression given in (|5.3p except the 
phase e^"^. Following [57|, using the famous formula 

det(l + X) = exp{Tr ln{l + X)} = exp |^(-l)'-iTr 

and the Schur polynomials given in Appendix[Dl we have an expression for det(l+ptp) in (IXT^ 
as 



det(l = Y^Sdx), XI ^ (-l)'-^Tr[(pM']//, 

1=0 

Xi = Si{x)=0 (/>M + 1). (J.l) 
Using the formula 

[det(l + X)]-5 =exp{Tr ln(l + X)-^} = exp |^(-1)' Tr (X')/(2/)| , 

the vacuum function ^M,Mip,P* ,t) (|3.12p is also expressed in terms of the Schur polynomials as 

= Tr[(ptp)]7(2/), 

(J.2) 



1=0 

^i = Si{0 = (/>M + 1). 



Rowe et al. have showed the number-projected SO{2N) wave function satisfies some recursion 
relations. They express it with the aid of the relations in a form of determinant [88] which is well 
known as the completely anti-symmetric Schur function in the theory of group characters [89] . 
In the present U{N) case, equation (IJ.ip is also given by a determinant form 



= ^det 



XI 1 
2X2 Xi 
3X3 2x2 




2 

Xi 





3 



IXi il-l)xi-i {l-2)xi-2 il-3)xi-3 



I - 1 
Xi 



i-iysii-x), 



which is exactly the same form as that given in [86]. The Schur function (fiix) satisfies the 
recursion relation and the differential formula 



^lix) 



l-i 



d 
dxi 



■^z(x) = (-l/+V«-^'(x). (J.3) 



By using the second equation of ()J.3p . we can rewrite the above recursion relation as 
1 ' 

^lix) = y Y.^-lY+^l'xi'^i-1'ix) {V>0 = 1). 



v=i 



With the aid of the second equation in p.ip with p = tan | • Iat and the explicit form of the 
Schur polynomials given in Appendix [Dl the xi Si{x) take simple forms, respectively, as 



X. = (-l)-)(tan^|)'-iV, ^.(X) = ,(^_^), 



1^-' (tan^l)'. 



(J.4) 
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Then using the first equation in (jj.ip and ()J.4p . we have 

oo N J. J. 

det(l +ptp) = = E TTT^rnTl'""' f )' = + f = D'"'' ' 

«=o z=o 

which coincides with the result by the direct calculation of det(l + p^p). The vacuum func- 
tion (IJ.2P is obtained as 

^m,m(p,P*, r) = (cos l)''e-'^\ eai^M,MiP,P* ,r) = {^M and ^M)- 
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